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Abstract 

Accurate supersymmetric spectra are required to confront data from direct and indirect searches of super- 
symmetry. SuSeFLA\Qis a numerical tool which is capable of computing supersymmetric spectra precisely 
for various supersymmetric breaking scenarios applicable even in the presence of flavor violation. The pro- 
gram solves MSSM RGEs with complete 3x3 flavor mixing at 2-loop level and one loop finite threshold 
corrections to all MSSM parameters by incorporating radiative electroweak symmetry breaking conditions. 
The program also incorporates the Type-I seesaw mechanism with three massive right handed neutrinos at 
user defined mass scales and mixing. It also computes branching ratios of flavor violating processes such 
as Ij — > Uj, Ij — > 3 ^i, 6 — > S7 and supersymmetric contributions to flavor conserving quantities such as 
(5^ — 2). A large choice of executables suitable for various operations of the program are provided. 
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at a high scales specified by the mechanism of supersymmetry breaking. RG equations including full inter- 
generational mixing are then used to evolve these parameters up to the electroweak breaking scale. The low 
energy supersymmetric spectrum is calculated at the scale where successful radiative electro-weak symmetry 
breaking occurs. At weak scale Standard Model fermion masses, gauge couplings are determined ncluding 
the supersymmetric radiative corrections. Once the spectrum is computed, the program proceeds to various 
lepton flavor violating observables (e.g., BR{^ — > 67), BR{t — > /i7) etc.) at the weak scale. 

Solution method: Two loop RGBs with full 3x3 flavor mixing for all supersymmetry breaking parame- 
ters are used to compute the low energy supersymmetric mass spectrum. Adaptive step size Runge-Kutta 
method is used to solve the RGBs numerically between the high scale and the electroweak breaking scale. 
Rerative procedure is employed to get the consistent radiative electroweak symmetry breaking condition. The 
masses of the supersymmetric particles are computed at 1-loop order. The third generation SM particles 
and the gauge couplings are evaluated at the 1-loop order including supersymmetric corrections. A further 
iteration of the full program is employed such that the SM masses and couplings are consistent with the 
supersymmetric particle spectrum. 

Restrictions: 

Unusual features: 

Additional comments: 

Several executables are presented for the user. 

Running time: 0.2 seconds on a Intel(R) Core(TM) i5 CPU 650 with 3.20 GHz. 



1. Introduction 

Low energy supersymmetry [T] is currently being probed by the Large Hadron Collider (LHC) at CERN 
and the Tevatron collider at the Fermilalj^ On the other hand, there is already a huge amount of informa- 
tion which has been collected and is being collected which gives information on low energy supersymmetric 
Lagrangian indirectly. For example, the flavor experiments in the hadronic and the leptonic sectors place 
strong constraints on the flavor off-diagonal entries in the Lagrangian. Similarly, the astrophysical data 
on dark matter which has been improved with the latest WMAP 7-year results [2] also strongly restricts 
the parameter space where the mass and the couplings of the lightest supersymmetric particle (LSP) corre- 
spond to the observed relic density. While experimental evidence for supersymmetry is definitely far more 
superior compared to the indirect detection of supersymmetry, the power of indirect experimental data to 
constrain the parameter space cannot be underestimated. Furthermore, as it is well known even if there is a 
positive experimental signal at the LHC, it would be very hard to reconstruct the supersymmetric breaking 
Lagrangian unambiguously due to the large number of degeneracies present in the parameter space which 
can give similar signals at the colliders [3]. It has also been noted that flavor violating observables and dark 
matter could help to break these degeneracies [4]. 

Obviously, the flavor observables depend on the supersymmetric model in which they are calculated. 
Unfortunately, most of the present supersymmetric mass spectrum calculators do not take in to consideration 
the effect of flavor violation in the running of soft mass parameters either in the hadronic sector or the 
leptonic sectoij^ In hadronic sector, typically the CKM is considered to be the only source of flavor violation, 
while this works very well, as long as one restricts to the scheme of Minimal Flavor Violation (MFV), in 
a more general scheme of supersymmetric models, such an assumption cannot be supported. Recently 



^The Tevatron collider has recently stopped functioning. 
^The recent version of SPheno |38| is an exception. 
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'observed' deviations from the Standard Model CKM paradigm O H] might find explanations in terms 
of a supersymmetric standard model with some amount of flavor violation ^ . To study the associated 
phenomenology of such kinds of models either for dark matter relic density, collider searches or other 
theoretical aspects such as threshold corrections to fcrmion masses, gauge coupling unification etc., would 
require precise computation of the mass spectrum in these models. 

In the leptonic sector, the case for flavor violation is even more stronger. Firstly, neutrinos have non-zero 
masses and secondly their flavor mixing is large as has been observed in the neutrino oscillation experiments. 
Most of mechanisms of generating neutrino masses and mixing inevitably lead to large flavor violation in 
supersymmetric theories. One of the simplest ways to generate neutrino masses is the so called 'seesaw 
mechanism'. In the present work, we have restricted ourselves to Type-I seesaw mechanism, though the 
program can be generalized to incorporate other seesaw mechanisms by adding the corresponding RGBs. 
The Type-I seesaw mechanisnj^has three singlet right handed neutrinos added to the Standard Model which 
leads to two additional terms, the Dirac mass term combining the left and right neutrino fields and the lepton 
number violating Majorana mass term for the right handed singlet fields. The interplay between these two 
terms leads to small Majorana masses to the left handed neutrinos in the limit of heavy Majorana masses 
for the right neutrinos. The supersymmetric version of the seesaw mechanism was proposed long ago [S] 
and some of its consequences for leptonic flavor violations have immediately been noticed. Over the years, 
other theoretical/phenomenological consequences of having right handed neutrinos has been observed. In 
the following we list some of them. 

• Yj, — Unification: The presence of right handed neutrinos could significantly modify regions where 
T-b Yukawa couplings unify at the GUT scale in unified theories like SO(IO) or SU(5) [101 [HI [H US] • 
This is due to the fact that the neutrino Dirac Yukawa couplings enter the renormalization group 
equations (RGE) of the Yh and Yr at 2-loop level and Yt at the 1-loop. This is enough to change the 
Yfc/Fr ratio at the weak scale, if the neutrino Yukawa couplings are large. 

• Lepton Flavor Violation: As mentioned previously, one of the main consequences of the seesaw 
mechanism in supersymmetric theories is the violation of lepton numbers leading to rare flavor violating 
decays [9] . This flavor violation will be generated through the RGE even if the supersymmetry breaking 
mechanism at the high scale conserves flavor as in mSUGRA. In particular GUT models, the generated 
flavor violations could be large enough to strongly constrain observability of supersymmetry at the 
LHC [HITniinilTH]. On going experiments like MEG and future experiments hke PRISM/PRIME 
and Super-B factories have enhanced sensitivity to large amounts of parameter space even with small 
mixing and small tan /3 [161 EH] • 

• Dark Matter: One of the most surprising phenomenological aspects of seesaw mechanism and SUSY- 
GUT models has been the impact on Dark matter phenomenology. The presence of a single right 
handed neutrino with a large Yukawa coupling could significantly enhance the efficiency of the elec- 
troweak symmetry breaking and thus making the focus point region unviable within mSUGRA like 
models i20j . Similarly, GUT effects can significantly effect the stau co-annihilation region [TH] . 
The co-annihilation region and the focus point regions seem to be most vulnerable to these effects in 
other GUT models and mSUGRA incorporating Type-II or Type-Ill seesaw mechanisms pTj. It has 
been recently realized that even flavor effects could play a role in the relic density calculations in the 
early universe [22] . 

• Hadronic flavor violation: Incorporating Type-I seesaw mechanism in Grand Unified Theories 
(GUT) also has effects on the hadronic sector. For example, CP violation in the neutrino sector 
could be transmitted to the quark sector in SU(5) or SO(IO) theories [23l[24|. The large phases of 
the neutrino mass matrix can be transmitted to the hadronic sector with effects in K and B physics 
phenomenology. More generally in GUT theories, hadronic and leptonic flavor violations are related 
to each other by the underlying GUT symmetry [25] . 



^For a summary of the seesaw mechanisms, please see 8^. 
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• Collider Signals: Lepton flavor violation which might be inevitable due to the presence of a seesaw 
mechanism, can lead to flavor violation in the sleptonic sector as we have mentioned above. Such 
flavor violation can be studied at the colliders by measuring the mass differences between the sleptons 
by observing such as sleptonic oscillations etc. [251 [?fl [251 [^ . 

• Gauge Coupling Unification: Finally, let us note that it has been pointed that the presence 
of massive right handed neutrinos with large yukawa couplings in the MSSM Lagrangian improves 
the accurate unification of the gauge couplings because as gets contribution from the right handed 
neutrinos through the RGE running at the two loop level [50] . 

The above points provide enough justification to determine the supersymmetric mass spectrum in seesaw 
models at a high precision level including the effects due to flavor violatiorj^ . Unfortunately, while there exist 
very good spectrum calculators for supersymmetric theories like ISASUSY \3T , SuSpect [34] and SOFTSUSY 
[35j . they do not consider full flavor violating structure in the computation of the soft spectrum cither 
in their RGB's or their mass matrice^ For these reasons, present versions of these programs might not 
be suitable for attacking problems listed above unless one significantly modifies them. Our program was 
written to address this deficiency in publicly available codes. We, however point out that the recent version 
of SPheno |39j is very similar to the program we are presenting here. It has full flavor structure for the soft 
masses as well as Yukawas at the 2-loop level and considers the full 6x6 mass matrices for the sparticles. 
The preliminary version of our program was flrst presented at [40] . The present version is an expanded and 
more structured version of the same. This paper explains the code in detail. We attempted to link the 
physics discussion with the file structure of the code wherever possible such that the user can modify the 
code with minimal efforts. 

SuSeFLAV is a program written in FortrcLn 95 in a fixed length format. The recommended compiler for 
this program is gf ortran available in various distributions of Linux. The program can be made executable 
using other Fortran compilers too, such as ifort, by modifying the Makefile in the main directory. In 
addition to studying the spectrum of supersymmetric particles the code also computes leptonic flavor vio- 
lating decays and some hadronic decays. We have implemented the SLHA 2 . 2 ^Tl format for dealing with 
the input parameters and output data. This way, the output of the code can be fed in to other publicly 
available programs either for computing Dark Matter relic density or for calculating supersymmetric parti- 
cle decays and production cross-sections at LHC etc. While the main set of RGB's are written for Type-I 
seesaw mechanism, extending the program for other seesaw mechanisms or even other models would require 
coding the RGB's from respective models. However, other parts of the code, like mass spectrum, one-loop 
corrections etc., which are given in separate files in the source directory can be used to suit the user's needs. 

The rest of the paper is organized as follows. In section [2] we describe the MSSM Lagrangian with Type- 
I seesaw mechanism. In section [3j we discuss the various supersymmetry breaking scenarios considered in 
SuSeFLAV. In section[4jwe describe the calculation of the low energy supersymmetric spectrum implemented 
in the program. In section [5] we briefly describe about various theoretical constraints we impose on the 
sparticle spectrum in SuSeFLAV. In section [6] we discuss about the various low energy observables computed 
in the program. In section [7| we show the instructions on how to install and execute SuSeFLAV. In | Appendix] 
[A] we write down all the tree level mass matrices of the sparticles and in Appendix B| we outline all the 



one-loop radiative corrections to these parameters SuSeFLAV. In [Appendix C we graphically show the flle 
structure in SuSeFLAV. We close with a comparison chart of sparticle masses with other available codes in 
Appendix D 



^Programs like SuperIso l31l and SUSYFLAVOR I32I compute the hadronic flavor violating processes for the given super- 
symmetric spectrum at low energy with high precision. 

® However there are programs which consider flavor violation in their mass matrices. The program SPICE 36 computes 
mass matrices with full flavor violation without considering intergenerational mixing in the RGE. After this work has appeared 
we have been informed that SOFTSUSY also has a new version where flavor violation is considered \37\ . 
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2. Minimal Supersymmetric Seesaw Model 

Since 1998, ever increasing data from neutrino sector has firmly established that neutrinos have non-zero 
masses and that their flavors mix with two of their angles being close to maximal and the third angle is 
non-zero [42] . One of the elegant mechanisms to generate non-zero neutrino masses is through the seesaw 
mechanism [43], where right handed singlet fields are added to the Standard Model particle spectrum. These 
singlet neutrinos break lepton number typically at a scale much larger than the standard model scale through 
their Majorana masses. 

Supersymmetric version of the seesaw mechanism is straight forward extension of the Minimal Su- 
persymmetric Standard Model [T] by adding right handed neutrino superfields. The field content of the 
MSSMRN(MSSM -I- Right Handed Neutrinos) and their transformation properties under the gauge group 
^SM = SU(3)c (g> SU(2)l «) U(1)f is given as 

z.^: (1,1,0); 

d^:(3,l,iy, (1) 

where Q and L stand for the SU{2)l doublet quarks and leptons, u^, d'^, e"^ and v'^ stand for the SU{2)l 
singlet quarks and leptons, and the two Higgs doublet chiral superfields are denoted by and H^. At the 
seesaw scale and above, > Af|., the superpotential takes the form : 

+ Y'^^iy^L.Hu + ^iHuHd - ^MR^i^fi^!, (2) 

where i,j — {1,2,3} are generation indices. Note that the right handed neutrino Majorana mass matrix 
is diagonal. The SU(2) and SU(3) contractions are suppressed in the above Lagrangian. In the program, 
the seesaw scale is taken to be = ^R^- ^^^^ scale, the right handed neutrino mass matrix can be 
diagonalized by a rotation of the right handed neutrino fields. The Dirac Yukawa matrix V is defined at 
this scale in the basis where the right handed neutrino mass matrix is diagonal. SuSeFLAV considers the 
inputs at the scale Mn^ where right handed neutrino masses as well as the neutrino Dirac Yukawa coupling 
matrix. Finally it should be noted that in the present version of SuSeFLAV, the ordering of the right handed 
neutrino mass eigenvalues is taken as A/^j^ < Mji^ ^ Mji^. We have not included the option of inverted 
hierarchy for the right handed neutrinos in the present version. The program will abort if such a choice is 
made. 

Below the seesaw scale, once the right handed neutrinos are integrated out, we are left with the five 
dimensional operator defining the the light neutrino mass matrix as 

W - Y^d'iQ.Ha + YX.u'iQjHu + Yf^etL.H^ 

+ ^^i-HuHd + LiHuLjHu (3) 

The light neutrino mass matrix is given at the weak scale after the SU(2)^ x U(l)y breaking as 

M.^^ {Hlf (4) 

where A represents the right handed neutrino mass scale or the seesaw scale and is the vev of the iJ„ 
superfield. The five dimensional operator is renormalized from the seesaw scale to the electroweak scale. 
These corrections can be significant for inverse hierarchal and degenerate spectrum for neutrino masses 
[44) . The present version of the program does not contain the renormalization for the light neutrinos, we 



L 
Q 



l,2,--j; e=:(l,l,+l) 
3,2,+i); «^:(3,1,-^ 
1,2,+i); i?.:(l,2,-i 
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refer users to use the existing programs like REAP 44 . The main reason for not including the RG effects 
for light neutrinos has been the famous ambiguity in relating the light neutrino masses to the neutrino 
Yukawa couplings [45]. Furthermore, we are more interested in the effects of seesaw mechanism on the 
soft supersymmetric masses and couplings. However, we do provide the option for the users to define the 
Yukawa matrices in terms of the Casas-Ibarra 7?.-parameterization [45] . In this parametrization the neutrino 
Yukawa matrix Y"^ is defined in terms of light neutrino masses as well as the right handed neutrino masses. 
In 7?.-parametrization, Y"^ at the seesaw scale defined as 



where Upmns is the Pontecorvo-Maki-Nakagawa-Sakata matrix 46 and TZ is any 3x3 orthogonal matrix. 



The V^/jj-^ and V^rj^ are defined as 



diagonal (v/M^, (6) 
T^^m;, = diagonal (^^7, v^) (7) 




where ki, k,2 and K3 are light neutrino mass eigenvalues. It is important to note that TZ can be complex 
in nature also, but in the present version of SuSeFLAV we take TZ to be real orthogonal matrix. One can 
parametrize the TZ matrix in terms of 3 angles but in SuSeFLAV we have not parametrized the TZ matrix and 
left all the 9 elements of TZ as user defined input. Various other parametrization of the Y^ matrix have also 
been carried out in |47j . Users interested in using the R-parameterization can use sinputs-rpar . in in the 
main directory for their computations. 

In addition to the user defined Y" at the seesaw scale, we also provide two other choices for and 
Affl based on Grand Unified Models like SO(IO). Both these cases consider hierarchal masses for the light 
neutrinos {M.^)- These are 

1. CKM Case: 

In this case the Y'' and Mr are given as 

VcKM (8) 
ht) 

Diagonal [Afi^J ={Mr.„ Mr^, Mr^} = {lO", 10^ 10*^} GeV 

where hu,hc, ht are the Yukawa couplings of the up, charm and the top quarks and Vckm is the quark 
sector mixing matrix. 

2. PMNS Case: 

In this case the Y'' and Mr are given as 

/ hu \ 

Y'' = \ K \ UpMNS (9) 

V ht J 

Diagonal[Af^] = {Mr^, Mr^, MrJ = {lO", 10^, 10^} GeV 

where Upmns is the leptonic mixing matrix. Both Vckm and Upmns are defined in the SuSemain . f file 
in the src/ directory. Below the seesaw scale, < ^'^r^j the right handed neutrinos (and sncutrinos) 
decouple from the theory, and the model is defined by MSSM. 

The soft supersymmetric breaking terms in the seesaw enhanced MSSM are same as in the MSSM 
together with the additional terms involving the right handed sneutrinos. These include the mass terms for 
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the gauginos, mass squared terms for all the scalar particles and also the bilinear terms and trilinear terms: 

+ m^^ \Huf + m^^ \Hd\^ + m\^ifLj + mfc^i'^e^ + m^^^^^*z7<=j- + . . . 
+ B^^HuHd + BA/,,^fi^^ + h.c. 

+ ^Q^u'^^Ka + Af^Q.d'jHa + MjLi^'jHa + ^L^l^^H^ (10) 

We use the factorization A" = A^Y", for all the trilinear couplings at the weak scal^ As noted before, 
at the energies < M^^, the right handed sneutrinos also decouple from the theory, along with the right 
handed neutrinos. In the program, we decouple the right handed neutrinos sequentially at different scales : 
the heaviest right handed neutrino at AZ/jg and the second heaviest one at Mn^ and the lightest one atM^j . 

The right handed neutrinos can be easily removed from the model to recover MSSM without right handed 
neutrinos, by choosing V = 0. This automatically decouples the right handed neutrinos in the theory. An 
explicit option is also provided in the input files, where by turning on/off the parameter rhn, one can either 
include/remove right handed neutrinos in the theory. Finally, let us note that quantum effects above the 
scale of seesaw ^ ^^r^ ^^^^ make the mass matrix AIj^ non-diagonal. The running of the Majorana mass 
matrix can also effect the B^f term in the soft potential, described below, at the 1-loop level. The can 
have implications for flavor physics and EDMs, if it is large through finite terms. These effects j48] are not 
computed in SuSeFLAV. 

The complete set of 2-loop RGE for a general superpotential and MSSM are presented in 0^. For the 
supersymmetric seesaw model, we use the RGE's from [50] , 



3. SUSY Breaking Mechanisms 

Supersymmetry is broken spontaneously in a hidden sector and is then communicated to the MSSM sector 
through the 'messenger sector'. The messengers could be gauge interactions or gravitational interactions. 
The result of this communication leads to soft supersymmetry breaking terms in the MSSM. While the 



form of the soft Lagrangian Eq. (10 1 is itself not dependent on the mediation mechanism, the physical 
quantities i.e., the masses, the couplings etc., are determined in terms of few 'fundamental' parameters 
depending on the mediation mechanism. Popular among such supersymmetry breaking schemes are (i) 
minimal Supergravity (mSUGRA) (ii) Gauge Mediated Supersymmetry Breaking (GMSB) (iii) Anomaly 
Mediated Supersymmetry Breaking (AMSB) (iv) Gaugino mediation (v) Moduli mediation etc. among a 
host of other possibilities [ST]. In addition there could be variations within each of the above schemes. In 
the SuSeFLAV, we have in-built (i) mSUGRA and some of its variations: (a) Non-Universal Higgs Models 
(NUHM), (b) Non-Universal Gaugino Models (NUGM) and (c) Complete Non-Universal Model (CNUM) 
and (ii) Gauge Mediated Supersymmetry Breaking (GMSB) models. The corresponding input files are given 
in the examples/ directory. An input file with completely non-universal soft parameters is also presented 
for supergravity mediation, where the users can define the boundary conditions of their choice. For other 
supersymmetric breaking models, the users can modify appropriately the slha.in file and run the program 
accordingly. In the section below we describe the two supersymmetric breaking scenarios considered in 
SuSeFLAV. 



3.1. mSUGRA and its variations 

Supersymmetry is broken spontaneously in a hidden sector and the communicated to the visible sector 
through the gravitational interactions. If the supergravity Kahler metric is canonical in matter fields, the 
soft terms resultant after integrating out the supergravity multiplet (while keeping the gravitino mass fixed). 



'^The RGEs are however, defined in terms of A^, which sets the format for inputs at the high scale. In mSUGRA, we have 
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are universal in nature |52j . The property universal refers to the flavor space i.e. all the soft terms take the 
same value irrespective of the flavor at the mediation scale. In SuSeFLAV we have considered the mediation 
scale to be Mqut- At this scale, all the soft terms are determined by four parameters and the sign of the 
/Lt-parameter. 

1. At Mqut the gaugino masses are universal to a value M1/2, i-e. 

Mi(Mgut) = M2(Mgut) = Af3(MGUT) = Ml/2 (11) 

2. The scalar and the Higgs masses are given by the parameter mg at A/qut- 

m|_(MG(7T) = m|^^ {Mgut) ^ [Mgut) = ra\XMGUT) = mf XMgut) = ml 1 

mji^ {Mgut) = {Mgut) = (12) 

3. The trilinear couplings are given by the parameter Aq at A/qut 

^ {Mgut) = Mj {Mgut) = A^^- {Mgut) = 1., (13) 

To specify the spectrum at the weak scale, two more parameters need to be fixed. First is the ratio of the 
vacuum expectation values (vevs) of the two Higgs fields, tan/3 = Vu/vd- Second is a discrete parameter, 
the sign of /i or the Higgsino mass parameter. The magnitude of ^ is fixed by the radiative electroweak 
symmetry breaking mechanism which has been incorporated in the program. 

3.1.1. Non-Universal Models 

In models based on Grand Unified theories, it has been proposed that the strictly universal feature of the 
soft masses might not be valid and in fact some amount of non-universalities can enter in a model-dependent 
fashion. For example, in models where the hidden sector gauge kinetic function is no longer singlet under the 
GUT group, gaugino masses would become non-universal at the high scale. However the non-universalities 
enter in a predictive fashion when a particular gauge group is chosen, as the ratios of the gaugino masses are 
now fixed by the Clebsch-Gordan coefficients of respective decomposition. These ratios are well known for 
the GUT models based on various gauge groups, e.g., SU(5) [H |M1 IMl liS] , SO(IO) [571[SHl[5g. Without 
resorting to any particular model, we have incorporated non-universal gaugino mass scenario in SuSeFLAV, 
by considering 

Mi{Mgut) ^ M2{Mgut) ^ M3{Mgut)- (14) 

The user has the freedom of choosing any ratios among these three parameters at GUT scale. The corre- 
sponding input file is s input s-nugm. in. A second class of non-universality which has been incorporated 
in SuSeFLAV is for the Higgs [6^. It has been argued that since in Grand Unified theories (GUTs) like 
SO(IO), all the matter sits in a single representation where as the Higgs sits in separate representation, the 
universality of the soft masses need not include Higgs, especially when supersymmetry breaking mediation 
happens close to the GUT scale. It has also been realized that introducing such non-universality makes the 
fi parameter free and thus leading to completely different phenomenology at the weak scale especially for 
dark matter. Thus the boundary conditions at the high scale in our notation are given by 

ui'qXMgut) = uiI^^Mgut) = ni}^^{MGUT) = m|XA^G[/T) = m? (A/gc/t) = ml 1 

'"^hA^'^gut) = mlo ; mH^{MGUT) = ml^ (15) 

Note that we intermittently use the notation miQ and m2o for the Higgs mass parameters as defined above 
in the non-universal Higgs mass model. The input file for this case sinputs-nuhm. in. It should be noted 
that a negative sign for any soft-mass as input would mean a negative sign for that soft-mass squared in the 
program. 

In addition to these non-universal input files, a completely generic input file called sinputs-cnum. in 
is provided where the user can provide all the supersymmetric breaking parameters at the high scale 
(mSUGRA) without any assumptions on their structure. This input file is more suited for models with 
fiavor structure at the high scale. 
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3.2. GMSB 



The second class of supersymmetric breaking models incorporated in SuSeFLAV is Gauge Mediated Su- 
persymmetric Breaking (GMSB). As before, supersymmetry is broken spontaneously in the hidden sector, 
but now communicated to the MSSM sector through gauge interactions. The minimal set of models under 
this category goes under the name, Minimal Messenger Model (MMM) [6T] . In this model, a set of messen- 
ger superfields transforming as complete representations of SU(5) (d Gsm) gauge group and couple directly 
to a singlet field which parameterized the supersymmetry breaking in the hidden sector. Supersymmetry 
breaking is then transmitted to the MSSM through SM gauge interactions. The following superpotential 
represents the messenger sector coupling to the hidden sector 

W = XX<^,^i (16) 

where $ and $ are messenger sector superfields transforming as 5 and 5 of SU(5) and i runs for the number 
of messenger sector superfields, typically i = [1,5]. The X superfield representing the hidden sector is 
parameterized by the vacuum expectation values for both its scalar component (X) as well as its auxiliary 
component (Fx). Gauge interactions with the messenger fields lead to gauginos attaining masses at 1-loop 
which at the Messenger scale are given by 

(M_.) = ^^^i^^Agf-£-Wn4^); (17) 

where the Messenger scale Mmess ~ ^{^) and A — {Fx) / {X) . Here na{i) is the Dynkin index for the 
messenger pair $, $ and the sum runs over all the messengers in each group. 
The scalars attain their masses from the 2-loop diagrams. These are given as 

m^(M_) = 2A2 / Y.^a{€) ) y (18) 

\ mess J ^ \ / 

Where Ca is the quadratic Casimir invariant of the MSSM fields and the function g{x) and j(x) are defined 
in [6 2) and [63j respectively. The leading order contribution to the tri-linear couplings comes from the 2-loop 
diagrams but they are suppressed by an extra Q;/47r factor compared to the gaugino masses. Thus to an 
very good approximation we can take 

ktf^Muess)^^ (19) 

Thus in minimal GMSB model considered in SuSeFLAV we have the following 5 parameters as inputs 

tan/3, sign(^), Mmess-, A and n. (20) 

The input file s input -gmsb . in and mgmsbl . 1 . in and so on, are built-in files to specify the inputs of the 
GMSB model in the examples/ directory of SuSeFLAV. 

4. Calculation of Supersymmetric Spectrum 

Once the user chooses the particular model of supersymmetry breaking by typing in the various pa- 
rameters in the relevant input file, the program computes the spectrum at the weak scale, checks for the 
various direct and indirect search limits and computes the various observables including the flavor violating 
ones like — ?■ e -I- 7. The computation of the spectrum involves several complicated intermediate steps 
which has been already explained in detail by various existing programs [BH [65] . In SuSeFLAV we follow 
a similar approach in computing the spectrum, however including flavor mixing as well as couplings with 
right handed neutrinos. In Fig. ([l]), we have shown the flowchart of the computation of the spectrum and 
the observables in the program. We can summarize the computation in terms of three steps which however, 
are not independent of each other as the procedure involves significant number of iterations. 
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Integrate RGEs to AIr^, 
heaviest RHN mass scale, 
Define 



SM inputs at Mz scale 
F", y^, y% aem, a,. First 
Iteration: no threshold 
corrections to inputs 



Run to high scale, impose 
SUSY breaking Boundary 
Conditions 



Calculate Low Energy 
Observables 



Yes 
— «— 



No 



No 



Global Convergence: 
Does ^ at msuay^ J^lt^^h 
and JJt Converge ? 



Integrate down to each RHIN 
mass scale , and decouple 
y below corresponding Mr, 



Run to EWSB scale 



Check |/t[^ > with 1-loop 
effective potential 



Check for consistent EWSB, 
Does /i converge? 



"Yes 



Figure 1: Flowchart of the program 



• RGE evolution: Using the MSSM RGE for the Yukawa and gauge coupHngs, run aU the known SM 
parameters hke gauge couphngs, Yukawa couphngs up to the scale of supersymmetry breaking. In the 
case of mSUGRA, run to the scale where the gauge couplings corresponding to SU(2)i (52) and U(1)y 

^gi = \l\ 9y^ meet, this determines the GUT scale (A/gut)- At the SUSY breaking scale (Mqut in 
case of mSUGRA and Mmess in case of GMSB), with the user defined input parameters and using the 
GUT scale SM parameters as the boundary conditions, run all the MSSM RGEs including those for 
the soft terms all the way down to MgusY- For the initial run A/susy is defined to be 1 TeV. If seesaw 
mechanism is switched on, the program takes in to consideration the running of user defined neutrino 
Yukawa couplings between the seesaw scale and the supersymmetric breaking scale in both directions. 

• Radiative Electroweak symmetry breaking: The resultant soft parameters at the Msusy are used to 
check if they satisfy the tree level electroweak symmetry breaking conditions and compute the ^ 
parameter. The full one loop effective potential corrections are then computed using the 'tree level' fi 
parameter, which is then used to derive the 1-loop corrections to the parameter. This is repeated 
iteratively until the n parameter converges. 

• Convergence of the Spectrum: In the final step, we run all the soft terms to the scale Mz where 
corrections to the SM parameters are added. We compute the supersymmetric corrections to the SM 
gauge couplings and the third generation fermion masses (i, b and r). The resultant masses are fed 
in to the RGE routine as shown in Fig. ([I]) and the soft spectrum is evaluated and run to MZ scale. 
This full iteration is continued until the SM third generation fermion masses converges to user defined 
precision (usually 0(10"'^ — 10"''')). Once this masses get converged we calculate various low energy 
observables (e.g. BR(/i — > ej), BR(6 — >■ 57) etc.). In section [6] we have discussed the low energy 
observables SuSeFLAV calculates. 

In the following we describe each of these steps in more detail. 
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RGE Evolution 

The standard model fermion masses and gauge couplings are the inputs to the program at the weak 
scale, taken to be equal to the Z-boson mass Mz in the program. The parameters are divided in to 
two subsets depending on whether radiative corrections are added or not. The parameters, masses of the 
first two generations quarks and leptons, for which we do not add radiative corrections are put in the 
file src/stdinputs .h. The values for these masses are taken from PDG 2012 Most of the other 

parameters such as the leptonic mixing matrix, UpmnSj the hadronic mixing matrix (Vckm) are defined in 
the file src/SuSemain . f . The MS values of Z-boson mass, the pole masses of the top quark m^°''^{mt) and 
tau lepton, m^°^'^ [rrir) and the MS mass of the bottom quark mass, m^^(mb), are left as user defined inputs. 
The MS gauge couplings, electromagnetic, aem{Mz) and the strong coupling, as{Mz) are also considered 
as inputs and are contained in the input files. The MS inputs are converted to DR as the RGEs are written 
in the DR scheme. The conversion for the gauge couplings is given by 



a. 



DR 



{Mz) 



1 



a 



MS 



{Mz) 



1 



„DR 



{Mz) 



1 



{Mz 



1 

47r 



(21) 



The so defined is in turn used to define the DR values of the ai,2- 



ai{Mz) 



9l 



^^^^iMz) 



a2{Mz 



9l 



a 



DR 



{Mz 



Itt 3 cos2 Oy/ 4tt sin 9w 

In a similar fashion, the bottom mass is converted from the MS to DR using [67l [68] 

as{Mz) 23al{Mz) , 3a2(Mz) , 13ai(Mz) 



m°«(Afz)=™r(m,) 



MS/ 



1 



3tt 



727r2 



327r 



2887r 



(22) 



(23) 



where the coupling constants appearing in the parenthesis are their DR values. The masses of the tau lepton 
and top quark are converted from their pole masses using the following relations: 



m^^{Mz) = m?°'« 



1 3 f . MMz) 
1 - - I ai(Mz) 



^?^{Mz) = mr-A,?,f^ 



where A^^^ is given by 



37r 



247r 



87r2 



with A*, = 2 In 



{mr'VM, 



and 



as{mt) 



a^'^{Mz) 



1 + 3^a?R(Mz)A, 



(24) 



(25) 



(26) 



The DR corrected masses are used to define the 3x3 Yukawa matrices at the Mz scale which form the 
inputs to the RGE. 



V sin /3 



Diag[m„,TOc,TOt] • Vckm 



Y" = 



V2 



V cos /? 



wcos/3 
Diag[me, m^, m.^] 



(27) 



We use the above defined Yukawas to run the full 2-loop RGEs from the weak scale {Mz) up to the scale 
at which the two gauge couplings {gi and 172) unify with an accuracy of 1%. For GMSB scenario this scale 
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is set by the user as messenger scale or Afmcss- In the case of right handed neutrinos three intermediate 
scales get introduced in the theory. As mentioned in the section [2] we consider the seesaw scale to be the 
mass of the heaviest right handed neutrino, M/jg. At this scale we set the neutrino Yukawa (Y^) and run 
the RGEs with this Yukawa up to the high scale. At the high scale, depending on the model, we set the 
SUSY breaking boundary conditions and then run the RGEs down to the heaviest right handed neutrino 
mass scale i.e., M^^. Below this mass scale we decouple the heaviest right handed neutrino by setting its 
couplings in the neutrino Yukawa matrix to zero. From M^j, we run down the RGEs to the next heaviest 
right handed neutrinos i.e., M^^ and then form M^^ to M^^. Below each of these scale, i.e., Mr^ and 
M/jj^ , we decouple the corresponding right handed neutrino by setting their couplings to zero. From Mp^-^ 
we run the RGEs down to the scale AfgusY- For the first iteration we take a guess value for Afgusy which 
is 1 TeV. From the second iteration onwards the Electro- Weak Symmetry Breaking (EWSB) scale is set to 
the geometric mean of the two stop masses or ^rrij^ • . At this scale we check for the EWSB condition 
and then we calculate the supersymmetric spectrum. A schematic picture of the integration of the RGE's 
in mSUGRA with seesaw mechanism is summarized in Fig. In GMSB models, the high scale is the 

messenger scale, Afmcss instead of A/qut and the procedure of the integration is very similar. The seesaw 
mechanism can be incorporated in this class of models as long as heaviest right handed neutrino is lighter 
than the messenger scale (Af^g < Afmoss)- 



GUT 



SUSY Breaking Boundary Condition: mSUGRA, NUHM 



Mr, 



Mz 



Define 



Decouple Yy 
row-by-row 



GUT 



Mr, 



M 



R2 



Mr, 



Mgusy 

Calculate Spectrum 



M, 



Inputs: F",y'*,y% 
aem,as, tan /3 



1-loop SUSY corrections to inputs 



Figure 2: Pictorial representation of evolution of RGEs 



.2. Radiative Electroweak Symmetry Breaking 
The tree level EWSB conditions at the AfgusY scale are defined as below 

ImI^ = ^ [tan 2(i{m\^ tan /? - m]j^ cot /3) - Af|] 



B, = ^[-|,„+n4,+2|Mp] (28) 
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where rnj^ and are the RGE output at Msusy- For consistent electroweak symmetry breaking we 

require > (]^ The tree level and RGE output of the other SUSY soft masses is used to calculate 
the tree level spectrum as described in [Appendix A| Radiative corrections can however significantly modify 
the tree level value of fi. Using tree level sparticle spectrum, we calculate the radiative corrections to the 
higgs potential up to one-loop order as given by BPMZ [7D]. The tadpoles modify the rnj^ and m'jj^ as 

rriH^ m]j - — ; m]j mj^ ~ —. (29) 



With these radiatively corrected higges, using Eq. (28), we calculate the radiatively corrected . This is 



repeated iteratively until the convergence of reaches the desired accuracy (default value is 0(10" )). 
This accuracy level can be changed by changing the parameter tol in file src/ewsbiterate . f . There could 
be regions where the ^ parameter does not converge within a small number of iterations. In such regions, 
the program considers the parameter point as |/i|-non convergent. Once \fi\ has converged we calculate 
as 

B, = '^[ml^+ml^ + 2\^^\'] (30) 

It is important to note here that at AfsuSY, SuSeFLAV checks for D-flat directions in the potential as well as 
whether the potential is unbounded from below. It also checks for charge and color breaking minima. More 
details about these checks are discussed in the next section. Even if these conditions are not satisfied the 
program still proceeds to compute the spectrum however, a flag is raised and written in the output file. 

Once fjL is converged, the program uses it to compute complete one loop corrections to the sparticle 
spectrum. We follow the work of BPMZ [70] in computing these correction^ Corrections to neutralinos 
and charginos are evaluated at external momenta equal to Msusy and the corrections to sfermion masses are 
calculated at an external momenta equal to their pole mass as prescribed by BPMZ. In [Appendix B| we have 
discussed more about these threshold corrections to the sparticles. As mentioned before all the parameters of 
the code are considered real, including the diagonalizing matrices. In determining the neutral higgs masses 
the user has a choice to employ approximations for two loop which are mostly top-stop enhanced [7T] or full 
one loop tadpole corrections as described in BPMZ or full one loop together with leading order two loop 
corrections. We have implemented the two loop corrections due to Slavich et. al [75] in our version 1.2. 
This is the default choice for the Higgs spectrum in the present version. 

4-3. Convergence of the Spectrum 

In final step, the program evaluates the full one loop fiavor conserving supersymmetric threshold correc- 
tions to SM parameters. The parameters which are corrected are mj, Wf,, mr^ ctg, Oiem and sir? Ow- One 
loop corrected running masses are given by the following, 

A„,(Mz) = Ef + AQf^, m,(Mz) = mf'^ [1 + A„,(Mz)] (31) 



mbiMzlMSSM - -I : ABPMZ ^•^^) 

' nib 



mr{Mz) = m?°'^ [1 + j:BPMZ^ ^33^ 

The quantity S^^*^-^, one loop correction to the top quark mass is evaluated at external momenta equal 
to m^^{Mz)- Whereas, /S.^'^^^ and S^-P^^^ are evaluated in the limit of external momenta tending to 



^Further, there should not be any tachyons in the physical Higgs spectrum. 

^Current version of the program does not include flavor violating contributions from sleptons to all the 1-loop corrections. 
The sleptonic contributions are neglected in this version in the presence of flavor violation. 
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zero. The expressions of these parameters are described in the Appendix B The three gauge couphngs get 
corrected as below 

ai{Mz) = r\ 2n (34) 

3 (1 - Aq_^^) cos"^ 9w 

(1 - Aq^^) sm 9w 
asiMz) = ^ _ ^ (36) 
where A.a,^ and A^^ are one loop corrections to the electromagnetic and strong coupling described in the 



Appendix B Note that sin^ 9^ used in the above expressions is also radiatively corrected. Iterative method 
is implemented to correctly evaluate SUSY contributions to sin^ 9w- The above corrected ai{Mz), a2{Mz) 
and a^{Mz) and also the third generation SM fermions (mf, m^, TOt-) masses are used as the input for the 
next long iteration. The iteration continues until the SM third generation fermions, namely top, bottom 
and tau mass are converged to user defined precision (usually 0(10"^)). This precision can be changed 
by the parameter named spectrum tolerance defined in the input files. Once the SM fermion masses get 
converged the program proceeds to calculate the various low energy observables. Finally, lets note that if 
both LL and RR type leptonic fiavor violation is present, it could lead to corrections to lepton self energies 
|73) , which are not included in the present version of the code. 



5. Theoretical and Phenomenological Constraints 

The requirement of consistent evaluation of supersymmetric spectrum involves a check for theoretical 
constraints such as charge and color breaking minima (CCB), scalar potential unbounded from below (UFB) 
and efficient electroweak symmetry breaking at EWSB scale. 

With every iterative step of the program SuSeFLAV checks for CCB and UFB conditions at the tree level 



|74) . These conditions are governed by equations 37 38 and are simultaneously implemented while computing 



the tree level ^ parameter (checking for efficient EWSB, see section 4.2 for the complete description) 



CCB: i(ml^+m) + \^,\^ + ml^>\Af\^ (37) 
UFB : m]i^ + m]j^ + 2\^i\^ > 2 at g > Mewsb (38) 

It is important to note that the complete MSSM spectrum is still calculated even if CCB and UFB 
conditions are not satisfied. However, if the supersymmetric scalar potential has a charge and color breaking 
minima which is lower than electroweak minimum a warning flag CCB is generated. Similarly, if the complete 
supersymmetric scalar is unbounded from below a warning flag UFB is generated. 

Apart form the above described theoretical constraints the program also imposes additional phenomeno- 
logical constraints on the obtained spectrum. From the direct non-observation of charged dark matter in the 
Universe we require nif > m^o or LSP being neutral. Regions for which this condition is not true is excluded 
as f LSP regions. Consequently the flag is marked as LSPSTAU. Also, we require the the spectrum to be 
non-tachyonic. If tachyonic spectrum is encountered it is flagged as TACSPEC. The program also indicates 
the sector where the tachyon occurs, for example: a tachyon in sleptonic sector is marked as TACSLP. Some 
lower bounds on various sparticle masses that result from direct searches at colliders, e.g., the lightest Higgs 
mass mfi > 114.5 GeV, and the Chargino mass m^± > 103.5 GeV |75j are also incorporated in the program. 
Points failing to satisfy these bounds are flagged as LEPH and LEPC respectively. The present version does 
not include recent direct search limits from LHC [751. 
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6. Low Energy Observables 



Once the complete supersymmetric spectrum is obtained we compute the foUowing low energy observables 
in SuSeFLAV. 

• Fine Tuning: In MSSM the standard model masses and gauge couplings, e.g. Mz, mt, mi, etc. are 
function of the input parameters of the model, i.e., for mSUGRA mp, M1/2, Aq, tan/3 and sign(/i). 
Now, we can recast the EWSB Eq. (pSl) as below: 



Ml 



tan 2/3 {mjf^^ tan /3 — cot /?) 



(39) 



where, all the parameters on the right side of Eq. ( 39 ) are at the weak scale. Given that Mz is know 
at a few percent level, it can be seen that some amount of 'tuning' between the parameters in the right 
side is needed. Various measures [771 IZHl [ZHl HO] have been proposed to quantify the 'tuning' needed 
in some input parameter. In SuSeFLAV, we have followed Barbieri and Giudice [771 [HO], in evaluating 
the fine tuning for a given parameter Xi of the model as 



^-/(M^,\) — 



(40) 



Where Barbieri-Giudice function or f{M'^,Xi) is defined as 

/(M|,A,) = ^^ (41) 
We derive the fine tuning in Af| with respect to /i^ for which /{M^, fj?) takes the following form 



M| 



(tan^ 13 + 1) 
(tan2 /3 - 1)3 



Itan^ 



The fine tuning in nit with respect to /i^ or f(Ml, mt) is expressed as follows 



1 TO- 



(42) 



(43) 



Ha 



• Electro-Weak Precision Measurement: The electroweak observables p parameter and sin^ 9w are 
defined in |44| and |45] 



P - 



Mw 



Mz cos 9 



w 



Mi 



sm e^^i-^ 



(44) 
(45) 



At the tree level p parameter is unity and sin^ 9w — 0.2286. A deviation from these values are observed 
at one-loop level stemming from SM and supersymmetric corrections to W and Z boson masses. A 
measure of this deviation to p parameter is given below 



rwo) 



n 



WW 
/\//-2 



(0) 



and p 



1 



(1 - Ap) 



(46) 
(47) 



The correction to sin^ 6w is described in Appendix B A precise measurement of these parameters by 
LEP, SLC and Tevetron imposes stringent constraints on Ap, requiring Ap < 10"'^ for a physically 
viable spectrum. 
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• BR(fe — > 57) Constraint: Another sector where the effect of SUSY particles can be seen is the 
radiative flavor changing decay of bottom quark, 6 — > 57 |81j . In the Standard model, this decay 
is mediated by loops containing up sector quarks and W bosons. In supersymmetric theories, if 
CKM is the only source of hadronic flavor violation, additional contributions to 5 — ?> 57 process come 
from chargino loops (x*), stop (t) squarks loops and charged Higgs bosons (H^) loops. As in the 
perturbation theory the SM and supersymmetric correction appear at the same order, the measurement 
of BR(5 — > 57) is a very powerful tool for constraining the SUSY parameter space. In SuSeFLAV we 
have followed Bartl et al. [82], which includes Standard model NLO as well as MSSM LO contributions. 
Like the other observables, this branching ratio is calculated for all the valid points and written in 
the SLHA file. More precise evaluations of this process are available publicly in the [3T1 [351 [351 US]- 
However, given that we are concentrating on evaluation of flavor violating observables in the presence 
of flavor violation in the soft sector, we would like upgrade our computation of this processes in a 
future version to include gluino contributions which could become important in the presence of squark 
flavor violating terms. Currently, we refer our users to couple the output of our program to one of 
the above existing codes depending on the accuracy required and input parameters they are interested 
in. For example Superlso [3T] computes the process at NNLO. It has an added advantage that it can 
compute other B-physics observables also. As of now our program does not have an interface with 
Superlso, but one can link to it using one of its programs (cmssm.c etc ) to our code. 

• Anomalous Magnetic Moment of Muon ((7^,-2): The muon anomalous magnetic moment {g^~2) 
has been precisely measured by Muon [g — 2) collaboration [84l[85]. Due to supersymmetric particles 
in the loop, (5^ — 2) gets non-negligible correction apart from the SM contribution. The experimentally 
measured value of anomalous magnetic moment of muon is: 

^exv = {Qp.-'^) ^ (11659208 ± 6) x lO^^". (48) 

Whereas, the difference in theoretical prediction by SM and experimental value, i.e aJ^^P — aj^*^ — 
(28.7 ± 8.0) X 10^^° [53]. Here the difference arises because of the fact of different estimates of 
hadronic vacuum polarization contribution. The contribution from the supersymmetric particles to 
{g^ — 2) is through (x" — jl/f/e) loop or through (x^ — v) loop. Supersymmetric parameter space can 
be severely constrained [57] as a supersymmetric explanation for the discrepancy with the experiment 
prefers a 'light' supersymmetric spectrum. Moreover, both these contributions are tan (3 enhanced, so 
large values of tan /? (> 30) are more severely constrained. We have taken the expression for one-loop 
contribution, due to supersymmetric particles, to the (g^, — 2) from Hisano et al. [H]. For a given 
set of input parameter SuSeFLAV calculates the (g^ — 2) and writes it into the SLHA flle. Two loop 
contributions are not added in the present version which could be important in the very large tan/3 
regime [55] . One of the main reasons for not including these contributions is the absence of two loop 
computations including flavor violation in the sleptonic sector. The present version includes flavor 
violating contributions at the 1-loop level which become important when the seesaw is switched on. 
We hope to generalize the two loop contributions including flavor violations in a future version. 

• Lepton Flavor Violating Decays: Subject to satisfying phenomenological constraints at AfsuSY, 
we evaluate decay rates of rare lepton flavor violating processes. The decay rates and branching ratios 
for the following processes are computed n — >■ 67, r — >■ 67, r — ^ /i7, ^ — >■ 67, e"'"e~e~, 

— )■ ji^ , — > e+e^e^ and /i — e conversion rate in the nuclei. The search for charged 

lepton flavor violating decays can play a pivotal role in studying and discovering new physics beyond 
the standard model at TeV scale and above. The amplitudes of these processes are also sensitive to 
the left-right mixing in the sleptonic mass matrices in addition to the intergenerational mixing [94]. 
In calculating the amplitudes we closely follow the notations and expressions provided in [lH [15]. 
The contribution to /i — e conversion rate in the nuclei is through penguin diagram where 7, Z is 
exchanged and through box diagrams containing — U ~ Qu.d loops or x~ ~ '^i ~ Qu,d loops. The 
current experimental upper bound on various LFV processes is tabulated in Table [l] 
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Table 1: Present Experimental Bounds on LFV Processes 



Process 


Experiment 


Present upper bound 


BR(/^ ^6 7) 


MEG [89] 


2.4 X 10"^^ 


BR(/i -> eee) 


SINDRUM [90] 


1.0 X 10^12 


CR(^ ein Au) 


SINDRUM-II [HI] 


7 X 10-13 


CR(^-^ einTi) 


SINDRUM-II [12] 


4.3 X 10-12 


BR(t ^ fij) 


BABAR/Belle [93] 


4.4 X 10-^ 


BR(t ^6 7) 


BABAR/Belle [23] 


3.3 X 10-^ 


BR(t fi fi fi) 


BABAR/Belle [35] 


2.0 X 10-8 


BR(t eee) 


BABAR/Belle 93 


2.6 X 10-8 



In the current version of SuSeFLAV we use the tree level masses and full 6x6 slepton mixing matrix to 
calculate LFV observables. The output is written in the BLOCK SuSeFLAVLOWENERGY in SLHA format. 

7. Executing SuSeFLAV 

Instruction to compile and install the package is provided in README file. SuSeFLAV uses the Diag 
package by T. Hanh to diagonalizc the mass matrice^^ SuSeFLAV package produces three executables when 
compiled, namely susef lav, susef lavslha and susef lavscEin in the bin/ sub-directory. To compute the 
spectrum for a single point the usage of executables susef lav and susef lavslha is recommended. The 
user must modify the corresponding input files sinputs.in and slha.in where all the input parameters 
are specified. To scan the parameter space the usage of the executable susef lavsceui is recommended. The 
corresponding input parameters are in the input file sinputs_scan. in. 

7.1. Sample Input/Output 

In this section we provide examples of input files which can be used to run the program. SuSeFLAV has 
two input/output modes. This can be broadly classified into SLHA I/O interface and non-SLHA or traditional 
SuSeFLAV I/O interface. The directory examples/ contains a variety of input files with which the user can 
run different models. 

7.1.1. SLHA Interface 

The main source file is runslha.f. To use the SLHA interface the user should execute susef lavslha 
by modifying the corresponding input file slha. in. Note that the user must rename the required SLHA file 
as slha.in to use that particular file as input. To run the program type the following command at the 
terminal 

. /susef lavslha 

This sample input /output corresponds to point 10.1.5 96 with Type-I see-saw mechanism with maximal 
mixing in mSUGRA/CMSSM scenario. 



^"Note that LAPACK 1951 library is a another option for an earlier version of SuSeFLAV for diagonalizing mass matrices. It is 
available on the website for downloading. We do not recommend the old version. 



17 



Block MODSEL 




TT tJtSJ-KSKi li 1U.\J KJiKS -L 






1 


1 




TT b u.y X d. 






6 


1 




it nnaT"V PI pnTTiT" V"i r»l r>Ti 

TT LLCLJ. J:^ J. J.CLV^J. V-L>-/-I-CLU-L \J±± 




Block 1 


3MINPUTS 




H C! -1- o n /-I o v J^ ]Vrr\/^^n "i titmtH e 
TT O UdULLclX U. riUUcX J.I1UU. U o 






1 


1 . 279340000E+02 




if ELXpna. V, IJ on rlODclX V.rlZi^ 






2 


1.166390000E-05 




# G_F6r']iii 






3 


1.172000000E-01 




TT cL-Lpild._o V,riZj^ Oil nOUcLL 






4 


9.118760000E+01 




if M7 (-nnl o"^ 

TT i IZj VpU-L" / 






5 


4.230000000E+00 




# rabCmb') SM MSbar 






6 


1.727000000E+02 










7 


1.770000000E+00 




if TTThan i "nnl pi j 

tT 111 LiCl.U.\kJV_/J_C/ 




Block MINPAR 




if TnTiiTl" "rarnf^l" PT'*^ 

Tr J.JJ.L/ LI U LJCLX CLillC U 






1 






TT lUU 






2 


7.000000000E+02 




if Tn1 9 
TT nil z 






3 


l.OOOOOOOOOE+01 




cLUb 






4 


l.OOOOOOOOOE+00 




# sign(inu) 






5 


O.OOOOOOOOOE+00 




# AO 




Block ; 


3USEFLAV 


it 


Algorithm specific inputs 




1 


1 


.OOOOOOOOOE-03 


Jf 

TT 


spectrum tolereince 




5 


2 


.OOOOOOOOOE+00 


# 


2-loop running (1 = Hoop; 2 = 21oop) 




3 


1 


.OOOOOOOOOE+00 


it 

TT 


print Control (0 = do not print output; 


1 = print output) 


2 


1 


.OOOOOOOOOE+00 


TT 


Right handed neutrino (1 = yes; = no) 




4 


1 


.OOOOOOOOOE+00 


it 

TT 


RHN mixing : 1 = ckm; 2 =PMNS; 3 = user 


defined 


7 


1 


.OOOOOOOOOE+06 


Jf 

TT 


MRl, Lightest rhn decoupling scale 




8 


1 


.OOOOOOOOOE+09 


it 

TT 


MR2, second lightest rhn decoupling scale 


9 


1 


.OOOOOOOOOE+14 


# 


MR3, Heaviest rhn decoupling scale 




10 





.OOOOOOOOOE+00 


# 


Dirac neutrino mixing matrix 1 , 1 




11 





.OOOOOOOOOE+00 


# 


Dirac neutrino mixing matrix 1,2 




12 





.OOOOOOOOOE+00 


# 


Dirac neutrino mixing matrix 1,3 




13 





.OOOOOOOOOE+00 


# 


Dirac neutrino mixing matrix 2,1 




14 





.OOOOOOOOOE+00 


# 


Dirac neutrino mixing matrix 2,2 




15 





.OOOOOOOOOE+00 


# 


Dirac neutrino mixing matrix 2,3 




16 





. OOOOOOOOOE+00 


# 


Dirac neutrino mixing matrix 3,1 




17 





.OOOOOOOOOE+00 


# 


Dirac neutrino mixing matrix 3,2 
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.OOOOOOOOOE+00 


# 


Dirac neutrino mixing matrix 3,3 





The corresponding SLHA output is generated in slha.out 



# Spectrum Output in SUSY Les Houches Accord 2 

# SuSeFLAV vl.2.0 

# D. Chowdhury, R. Garani and S. K. Vempati, hep-ph/1 109 . 3551 

# For bug reports or any other queries please send email to suseflav@cts.iisc.ernet.in 

# Created on 12.06.2012 at 17:30 Hrs 
# 

BLOCK SPINFO # Program information 

1 SuSeFLAV # Spectrum calculator 

2 1.2.0 # Version number 

# 

BLOCK MODSEL # MODEL NAME 
1 1 # mSUGRA 

3 4 # mSUGRA + RHN 

6 3 # Lepton and quark flavor is violated 
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BLOCK MINPAR # Input parameters 



1 


1, 


.75000000E+02 


# 


mO 


2 


7, 


.OOOOOOOOE+02 


# 


m_l/2 


3 


1. 


.OOOOOOOOE+01 


# 


tanbeta(mz) 


4 


1, 


,00000000E+00 


# 


sign(mu) 


5 


0, 


.OOOOOOOOE+00 


# 


AO 



# 

BLOCK 



SMIMPUTS 
1 
2 
3 
4 
5 
6 
7 



# Standard Model inputs 
1.27934000E+02 # alpha_em (M_Z)~MSbar 



1 . 16639000E-05 

1.17200000E-01 
9.11876000E+01 
4.23000000E+00 
1.72700000E+02 
1.77000000E+00 



# G_F [GeV--2] 

# alpha_S(M_Z)~MSbar 

# M_Z pole mass 

# mb(mb)"MSbar 

# mt pole mass 

# mtau pole mass 



BLOCK EXTPAR 



n 
u 


1, 


.97682880E+16 


# 


Unif i. 




# Mass Spectrum 






ru\a coae 




ni3.s s 


particli 


OA 


8, 


05095481E+01 




w+ 


OR 


± ■ 




# 


h 


oo 


9. 


.79136838E+02 


# 


H 


OD 


9, 


.78810446E+02 


# 


A 


QV 
O ( 


9. 


.82557324E+02 


# 


H+ 


1000001 


1, 


. 29049 183E+03 


# 


~d_l 


1000003 


1, 


,33317083E+03 


# 


~d_2 


1000005 


1. 


. 33782720E+03 


# 


~d_3 


2000001 


1, 


. 33783952E+03 


# 


~d_4 


2000003 


1, 


.39947428E+03 


# 


~d_5 


2000005 


1. 


.39948751E+03 


# 


~d_6 


1000002 


1, 


.08957331E+03 


# 


~u_l 


1000004 


1, 


,32515013E+03 


# 


~u_2 


1000006 


1. 


. 34374200E+03 


# 


~u_3 


2000002 


1, 


. 34410204E+03 


# 


~u_4 


2000004 


1, 


.39733996E+03 


# 


~u_5 


2000006 


1. 


. 39738676E+03 


# 


~u_6 


1000011 


3, 


. 12007210E+02 


# 


~1_1 


1000013 


3, 


, 18856813E+02 


# 


~1_2 


1000015 


3. 


.18880311E+02 


# 


-1_3 


2000011 


4, 


. 94924639E+02 


# 


~1_4 


2000013 


4, 


,96469065E+02 


# 


~1_5 


2000015 


4. 


.96470193E+02 


# 


-1_6 


1000012 


4, 


.87191917E+02 


# 


~nu_l 


1000014 


4, 


,90375557E+02 


# 


~nu_2 


1000016 


4, 


.90381909E+02 


# 


~nu_3 


1000021 


1, 


.57599127E+03 


# 


'g 


1000022 


2, 


,89580143E+02 


# 


~chi_ 


1000023 


5, 


.44321259E+02 


# 


"chij 


1000025 


8. 


. 58605742E+02 


# 


~cM_: 


1000035 


8. 


. 6943747 lE+02 


# 


~chi_' 



10 
20 
30 
40 
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1000024 5.33851586E+02 # ~chi_l+ 

1000037 8.62050764E+02 # ~cM_2+ 

# 

BLOCK NMIX # Neutralino Mixing Matrix 



# 



1 

J. 


1 


Q 




Jf 

TT 


M 


1 1 




9 


-•i 
J. 


n'^fiRPj^^RTT— no 

. UODOOOOOIli 


Jf 
TT 


1M_ 


1 


1 


3 


6 


. 20427374E-02 


# 


N. 


.13 


1 


4 


-2 


.65332008E-02 


# 


N. 


.14 


2 


1 


2 


.32098394E-02 


# 


N. 


.21 


2 


2 


9 


.80197789E-01 


# 


N. 


.22 


2 


3 


-1 


.61409861E-01 


# 


N. 


.23 


2 


4 


1 


.12340793E-01 


# 


N. 


.24 


3 


1 


2 


.46412419E-02 


# 


N. 


.31 


3 


2 


-3 


.55274833E-02 


# 


N. 


.32 


3 


3 


-7 


.04937074E-01 


# 


N. 


.33 


3 


4 


-7 


.07950795E-01 


# 


N. 


.34 


4 


1 


-5 


.92863909E-02 


# 


N. 


.41 


4 


2 


1 


.94531848E-01 


# 


N. 


.42 


4 


3 


6 


.87867194E-01 


# 


N. 


.43 


4 


4 


-6 


.96764815E-01 


# 


N. 


.44 


.OCK UMIX 


# Chargino Mixing 


Matrix U 


1 


1 


9 


.74989137E-01 


# 


U. 


.11 


1 


2 


-2 


.22252520E-01 


# 


U. 


.12 


2 


1 


2 


.22252520E-01 


# 


U. 


.21 


2 


2 


9 


.74989137E-01 


# 


U. 


.22 



# 



BLOCK VMIX # Chargino Mixing Matrix V 
1 1 9.88224902E-01 # V_ll 

1 2 -1.53008313E-01 # V_12 

2 1 1.53008313E-01 # V_21 
2 2 9.88224902E-01 # V_22 

# 

BLOCK USQMIX # squark Mixing Matrix 



1 


1 


1. 


. 10996359E- 


-03 


# 


USQMIX_11 


1 


2 


-5, 


.26653686E- 


■03 


# 


USQMIX_12 


1 


3 


3, 


,55285163E- 


■01 


# 


USQMIX_13 


1 


4 


8. 


. 74704256E- 


■04 


# 


USqMIX_14 


1 


5 


-4. 


.21826454E- 


■03 


# 


USqMIX_15 


1 


6 


9, 


.34732542E- 


■01 


# 


USQMIX_16 


2 


1 


1. 


.12595131E- 


■02 


# 


USQMIX_21 


2 


2 


-5, 


.57113248E- 


■02 


# 


USQMIX_22 


2 


3 


9, 


,23559656E- 


■01 


# 


USQMIX_23 


2 


4 


2. 


. 86784533E- 


■02 


# 


USQMIX_24 


2 


5 


-1, 


. 38092429E- 


■01 


# 


USQMIX_25 


2 


6 


-3. 


.52015712E- 


■01 


# 


USQMIX_26 


3 


1 


2. 


. 97503603E- 


■05 


# 


USC!MIX_31 


3 


2 


-9, 


.98459110E- 


■05 


# 


USQMIX_32 


3 


3 


1, 


.03612223E- 


■03 


# 


USQMIX_33 


3 


4 


9. 


. 77509497E- 


■01 


# 


USC!MIX_34 


3 


5 


2. 


. 10888527E- 


■01 


# 


USC!MIX_35 


3 


6 


-3. 


. 57456096E- 


■04 


# 


USC!MIX_36 





1 

± 


9 


1 H71 ^n^RF- 
. XO ( lOUOOCi 


uo 


# 


TjqnMTY 


_^x 


A 
*± 


9 


9 






TT 


TTonMTY 


AO 


A 
*± 


o 


J. . 






Jt 
TT 


TTCinMTY 
Uol^rlXA_ 




A 
•± 


•± 


_9 

z . 




W X 


# 

Tr 


TjqnMTY 


_*±'± 




c 
O 


Q 


. D ( Drt f ZOOJZi 


U X 


a 

Tr 


iTqnMTY 
uoi^rix A_ 


.^0 


A 
*± 


p. 
D 


. 




-09 
uz 


Jt 
Tr 


TjqnMTY 

Uol^rlXA_ 




C 

o 


1 


q 


. OZZXZQV/ZCi 


UX 


jf 
Tr 


TjqnMTY 


_ox 


c 

o 


z 


•1 

J. , 


R77R'^7nnF- 


-01 

\J X 


# 
Tr 


TjqnMTY 
uoi^nx A_ 


R9 
.oz 


5 


3 


-9, 


.34002684E- 


■04 


# 


USQMIX. 


.53 


5 


4 


5. 


.08408180E- 


■04 


# 


USQMIX. 


.54 


5 


5 


-2, 


401 OR*^! QF- 


■03 


# 


U Oi^l IX A_ 


RR 


5 


6 


2, 




■04 


# 


TTCIOMTY 


Rfi 
.00 


6 


1 


-1. 




■01 


# 




.0 J- 


6 


2 


9, 


RORI fi99QF- 

, OL/D XOj^Zyili 


■01 


# 


U 11 X A _ 


fi9 


6 


3 


5, 




■02 


# 




.DO 


6 


4 


2. 




■03 


# 


TTqnMTY 


fiA 


6 


5 


-9, 




■03 


# 


U 11 X A _ 


fiR 
.DO 


6 


6 


-1, 


. J. y woy / Hi 


■02 


# 


TTqOMTY 


.DD 


.OCK DSQMIX 


TT oULLdXA 


Mixing 


1 id UX X A 




1 


1 


7, 




■03 


# 


T^qOMTY 
XiOl^l IX A_ 


1 1 
.XX 


1 


2 


-3. 




■02 


# 


jjoi^rix A_ 


1 9 
_ xz 




3 


9 




■01 




L/Ol^riX A_ 


1 

. xo 


1 


4 


-1, 


nA'^7'^1 01 F- 


■05 


# 


L/oi^nx A_ 


1 A 
. x*± 


1 


5 


-2, 


4.1 fi79'^dRF- 


■04 


# 


nqnMTY 

iyOL^i IX A_ 


1 R 
. J- 


1 


6 


2. 


1 '^9R'^fi1 RF- 


■01 


# 


nqnMTY 
j-'Oi^rixA_ 


1 


2 


1 


-2, 


P1 AQn^7AF- 


■03 


# 


nqnMTY 
x/oi^nx A_ 


91 
_Z X 


2 


2 


1, 


'^RfiA'^QQ9F- 


■02 


# 


nqnMTY 

iyOl^nX A_ 


99 
.zz 


2 


3 


-2. 


1 9RfiRR'^0F- 


■01 


# 


nqnMTY 
jjoi^rix A_ 


9"^ 
_zo 


2 


4 


1, 


7Q977A1 ftF- 
. f yz ( ( fx OIL 


■05 


# 


nqnMTY 
x/oi^nx A_ 


9A 


2 


5 


9. 




■04 


# 


nqnMTY 

JJoL^l IX A_ 


9R 
.zo 


2 


6 


9. 


7fiQR979'^F- 


■01 


# 


TjqnMTY 

L/oi^nx A_ 


Of, 

.ZD 


3 


1 


6, 


1 0Q9fi9Q1 F- 


-06 


# 


nqnMTY 
jJou/rix A_ 


•^1 
.OX 


3 


2 


5, 




■03 


# 


nqnMTY 

L/oi^nx A_ 


•^9 
.oz 


3 


3 


6. 


PfiOAAfi99F- 


■04 


# 


r»qnMTY 

L/Ol^l IX A_ 


.00 


3 


4 


-2, 


RQ1 ^91 ^RF- 


■05 


# 


nqnMTY 

JJOlj^rlX A_ 


•^A 
.0*1 


3 


5 


9, 




■01 


# 


nqnMTY 

L/oL^nx A_ 


•^R 
.00 


3 


6 


-9. 


ARRRQQQ1 F- 


■04 


# 


nqnMTY 

JJOl^l IX A_ 


.OD 


4 


1 


2. 




■04 


# 


nqnMTY 

L/Ol^l IX A_ 


A1 

.tX 


4 


2 


-7, 


QRfiftAQfiAF- 


■07 


# 


nqnMTY 

XiOl^l IX A_ 


A9 

.tZ 


4 


3 


1. 


1 R9^9fiR'^F- 


■05 


# 


nqnMTY 

JJOl^i IX A_ 


A*^ 

."O 


4 


4 


9. 




■01 


# 


nqnMTY 

JJOl^l IX A_ 


AA 


4 


5 


2, 


RPQAA557fiF- 


■05 


# 


nqnMTY 


AR 


4 


6 


-1. 


. 50004569E- 


■05 


# 


DSQMIX. 


46 


5 


1 


-1. 


.75678721E- 


■01 


# 


DSQMIX. 


.51 


5 


2 


9, 


.83552600E- 


■01 


# 


DSQMIX. 


.52 


5 


3 


4. 


. 12376560E- 


■02 


# 


DSQMIX. 


.53 


5 


4 


5, 


.24615573E- 


■05 


# 


DSQMIX. 


.54 


5 


5 


-5, 


.82957991E- 


■03 


# 


DSQMIX. 


.55 


5 


6 


-5. 


. 17098306E- 


■03 


# 


DSQMIX. 


.56 


6 


1 


9. 


.84411220E- 


■01 


# 


DSQMIX. 


.61 


6 


2 


1. 


.75875037E- 


■01 


# 


DSQMIX. 


.62 



6 3 -1.15326259E-03 # DSC!MIX_63 

6 4 -2.91633802E-04 # DSqMIX_64 

6 5 -1.04177570E-03 # DSQMIX_65 

6 6 1. 44250 146E-04 # DSC!MIX_66 



BLOCK SELMIX # Slepton Mixing Matrix 



1 
1 


■1 
1 


1 . 


.05172917E- 


■05 




QTTT MTY 


1 1 


1 


O 

z 


— 


.07243990E- 


■05 


ff 


QPT MTY 
OIL i_irilA_ 


■1 9 


1 


O 


i 

1 , 


.03659288E- 


■01 




QPT MTY 


_ -LO 


± 


/I 
ft 


a, 
D . 


. 00795431E- 


■08 


ff 


QTTT MTY 


_-Lft 


1 


C 

O 


1 , 


.11136067E- 


■05 


4f 

ff 


QTT MTY 




1 


D 


Q 


.94612864E- 


■01 


ff 


QFT MTY 


1 

_ -LO 


o 


■1 
1 


O . 


.08963419E- 


■08 


ff 


QTTT MTY 


91 


o 
z 


Z 


C 

D , 


,24290695E- 


■03 


ff 

ff 


QFT MTY 
oCiLirll A_ 


99 
_ZZ 




o 
O 


_9 

z . 


.02732368E- 


■06 


ff 


QFT MTY 
OX-iijllX A_ 


_ZO 


n 
z. 


A 
ft 




.78386431E- 


■09 


ff 
ff 


QTTT MTY 


9A 


o 
z 


C 

O 


Q 


.99980513E- 


■01 


ff 

ff 


QFT MTY 
oJ_iI_irll A_ 


9C^ 
_ZO 




D 


i 

1 , 


.17032561E- 


■05 


ff 


cpT MTY 
o]_iJ_irll A_ 


_ZD 


O 


■1 
1 


9 
Z , 


.98889602E- 


■05 


ff 
ff 


QT?! MTY 




o 

O 


z 


D , 


.33346635E- 


■11 


ff 
ff 


cpT MTY 


_oz 


Q 
O 


O 


1 . 


.18627921E- 


■08 


ff 


QTTT MTY 
OEiLirll A_ 


_oo 


Q 
O 


ft 


i 

1 , 


.OOOOOOOOE+00 


ff 


c;T7T MTY 




O 
O 





C 
D , 


.78606456E- 


■09 


ff 


QFT MTY 
oJZiLirll A_ 


_oo 


Q 
O 


D 


D . 


.19572827E- 


■08 


4f 

ff 


CPT MTY 
oiliLiril A_ 


_oD 


/I 


■1 
1 


Q 


.84024127E- 


■03 


if 
ff 


CPT MTY 
OI-iLiril A_ 


A1 


ft 


9 
Z 


—A 

ft , 


.73545901E- 


■02 


4f 
ff 


c;T7T MTY 


A9 


4 


3 


9, 


. 93448538E- 


■01 


# 


SELMIX. 


.43 


4 


4 


-3, 


.12319765E- 


■07 


# 


SELMIX. 


.44 


4 


5 


2. 


. 98862230E- 


■04 


# 


SELMIX. 


.45 


4 


6 


-1, 


. 03540457E- 


■01 


# 


SELMIX. 


.46 


5 


1 


9, 


,99950259E- 


■01 


# 


SELMIX. 


.51 


5 


2 


2. 


.09161645E- 


■03 


# 


SELMIX. 


.52 


5 


3 


-9, 


.70065714E- 


■03 


# 


SELMIX. 


.53 


5 


4 


-2, 


,98872962E- 


■05 


# 


SELMIX. 


.54 


5 


5 


-1. 


.30585027E- 


■05 


# 


SELMIX. 


.55 


5 


6 


1, 


. 00054245E- 


■03 


# 


SELMIX. 


.56 


6 


1 


-1, 


.62739211E- 


■03 


# 


SELMIX. 


.61 


6 


2 


9. 


. 98856442E- 


■01 


# 


SELMIX. 


.62 


6 


3 


4, 


.71237978E- 


■02 


# 


SELMIX. 


.63 


6 


4 


4, 


.78864799E- 


■08 


# 


SELMIX. 


.64 


6 


5 


-6. 


. 23573698E- 


■03 


# 


SELMIX. 


.65 


6 


6 


-4, 


. 86038870E- 


■03 


# 


SELMIX. 


.66 


.OCK SNUMIX 


# Sneutrino 


Mixing Matrix 


1 


1 


4, 


.85609162E- 


■03 


# 


SNUMIX. 


.11 


1 


2 


-2, 


,34212181E- 


■02 


# 


SNUMIX. 


.12 


1 


3 


9. 


.99713892E- 


■01 


# 


SNUMIX. 


.13 


2 


1 


-5. 


.74301350E- 


■02 


# 


SNUMIX. 


.21 


2 


2 


9. 


.98069089E- 


■01 


# 


SNUMIX. 


.22 


2 


3 


2. 


.36616496E- 


■02 


# 


SNUMIX. 


.23 


3 


1 


9. 


.98337717E- 


■01 


# 


SNUMIX. 


.31 


3 


2 


5. 


.75286069E- 


■02 


# 


SNUMIX. 


.32 



3 3 



-3.50163123E-03 # SNUMIX_33 



# 

BLOCK ALPHA # Higgs mixing 

-1.05473355E-01 # 

# 

BLOCK HMIX Q= 1 . 20195104E+03 



Mixing angle in the neutral Higgs boson sector 



1 8.50892351E+02 

2 9.63860969E+00 

3 2.43818892E+02 

4 9 . 93075833E+05 



# DRbar Higgs Paraimeters 

# mu(Q) 

# tanbeta(Q) 

# vev(Q) 

# MA"2(Q) 



BLOCK GAUGE Q= 1 . 20195104E+03 # The gauge couplings 



1 3.62957251E-01 

2 6.42610014E-01 

3 1 . 04880243E+00 



# gprime(Cj) DRbar 

# g(Q) DRbar 

# g_3(Q) DRbar 



BLOCK TU q= 1.20195104E+03 # The trilinear couplings 



1 


1 


-1, 


.26688884E-01 


# 


TURG. 


.11 


1 


2 


-2. 


.93067837E-02 


# 


TURG. 


.12 


1 


3 


-3, 


.97143909E-04 


# 


TURG. 


.13 


2 


1 


2. 


.93036897E+01 


# 


TURG. 


.21 


2 


2 


-1, 


.26686864E+02 


# 


TURG. 


.22 


2 


3 


-5, 


.37734251E+00 


# 


TURG. 


.23 


3 


1 


-1, 


.11137936E+02 


# 


TURG. 


.31 


3 


2 


5, 


.35228445E+02 


# 


TURG. 


.32 


3 


3 


-1. 


.29451330E+04 


# 


TURG. 


.33 


.OCK TD 


Q= 


1.20195104E+03 


# The 


tr: 


1 


1 


-6. 


. 05492867E+00 


# 


TDRG. 


.11 


1 


2 


-1, 


.02056677E-04 


# 


TDRG. 


.12 


1 


3 


2, 


.38753240E-03 


# 


TDRG. 


.13 


2 


1 


-2. 


.04113446E-03 


# 


TDRG. 


.21 


2 


2 


-1, 


.21088577E+02 


# 


TDRG. 


.22 


2 


3 


-2, 


.29965104E-01 


# 


TDRG. 


.23 


3 


1 


1. 


.22702682E+00 


# 


TDRG. 


.31 


3 


2 


-5. 


. 90930668E+00 


# 


TDRG. 


.32 


3 


3 


-2. 


.95665819E+03 


# 


TDRG. 


.33 



BLOCK TE Q= 1 .20195104E+G3 

1 1 -1.44939639E-01 

1 2 -2.31067898E-08 

1 3 -4.08400006E-07 

2 1 -4.82631701E-06 
2 2 -3.02751406E+01 

2 3 4.10638455E-04 

3 1 -1.53360878E-03 
3 2 7.38111171E-03 
3 3 -5.16767721E+02 

# 

BLOCK YU Q= 1.20195104E+03 

1 1 6.44039888E-06 

1 2 1.48968552E-06 



# The trilinear couplings 

# TERG_11 

# TERG_12 

# TERG_13 

# TERG_21 

# TERG_22 

# TERG_23 

# TERG_31 

# TERG_32 

# TERG_33 

# The top Yukawa coupling 

# YU_11 

# YU_12 
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1 


3 


9 


.41182558E- 


■08 


# 


YU. 


.13 


2 


1 


_-| 


,48971566E- 


-03 


# 


YU. 


.21 


2 


2 


c 
D . 


.44041987E- 


-03 


# 


YU. 


.22 


2 


3 


9 


.72567759E- 


■04 


# 


YU. 


.23 


3 


1 


7 


,22768067E- 


-03 


# 


YU. 


.31 


3 


2 


_q 


.48078574E- 


-02 


# 


YU. 


.32 


3 


3 


o 
o . 


.43790071E- 


■01 


# 


YU. 


.33 


.OCK YD 


n— 


1.20195104E+03 


# The ( 


1 


1 


2. 


. 53400760E- 


■04 


# 


YD. 


.11 


1 


2 


-1, 


.10093360E- 


-09 


# 


YD. 


.12 


1 


3 


2, 


.56831364E- 


-08 


# 


YD. 


.13 


2 


1 


-2. 


.20186781E- 


■08 


# 


YD. 


.21 


2 


2 


5, 


,06812345E- 


-03 


# 


YD. 


.22 


2 


3 


-2, 


.47377375E- 


-06 


# 


YD. 


.23 


3 


1 


1. 


.31978305E- 


■05 


# 


YD. 


.31 


3 


2 


-6, 


.35600675E- 


-05 


# 


YD. 


.32 


3 


3 


1, 


.31856452E- 


-01 


# 


YD. 


.33 



,OCK YE 


Q= 


1.20195104E-I-03 


# The ■ 


1 


1 


2, 


.77703519E- 


-05 


# 


YE. 


.11 


1 


2 


-1. 


. 17975032E- 


-11 


# 


YE. 


.12 


1 


3 


4, 


.03483137E- 


-10 


# 


YE. 


.13 


2 


1 


-2, 


.46439380E- 


-09 


# 


YE. 


.21 


2 


2 


5, 


.80081939E- 


-03 


# 


YE. 


.22 


2 


3 


-4. 


. 05872732E- 


-07 


# 


YE. 


.23 


3 


1 


1, 


.45732796E- 


-06 


# 


YE. 


.31 


3 


2 


-7, 


.01783388E- 


-06 


# 


YE. 


.32 


3 


3 


9. 


. 96765043E- 


-02 


# 


YE. 


.33 



BLOCK MSOFT Q= 

1 

2 

3 
21 
22 
31 
32 
33 
34 
35 
36 
41 
42 
43 
44 
45 
46 
47 
48 
49 

# 



1.20195104E-I-03 
2.95504689E-H02 
5.45936065E+02 
1.53262491E-^03 
2.14879878E-I-05 
6 . 96254426E+05 
4 . 94404728E-I-02 
4. 943968 17E-I-02 
4.91242768E-H02 
3.15861187E-H02 
3.15846380E-H02 
3.11472426E-F02 
1.39829222E-H03 
1.39811849E-H03 
1.29136113E-H03 
1.34421791E-^03 
1 . 34420302E-H03 
1.11340022E-I-03 
1.33760077E-H03 
1 . 33759064E-H03 
1.33101944E-H03 



# soft SUSY breaking masses at scale Q 

# M_l 

# M_2 

# M_3 

# M-2_Hd 

# M-2_Hu 

# M_eL 

# M_muL 

# M_tauL 

# M_eR 

# M_muR 

# M_tauR 

# M_qlL 

# M_q2L 

# M_q3L 

# M_uR 

# M_cR 

# M_tR 

# M_dR 

# M_sR 

# M_bR 
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BLOCK MSQ2 Q= 1 . 20195104E+03 # n~2_Q 

1 1 1.95522114E+06 # mSQRG_ll 

1 2 1.03013845E+02 # mSC!RG_12 

1 3 -2.38957910E+03 # mSQRG_13 

2 1 1.03013845E+02 # mSQRG_21 
2 2 1.95473530E+06 # mSQRG_22 

2 3 1 . 15080756E+04 # mSQRG_23 

3 1 -2.38957910E+03 # mSqRG_31 
3 2 1.15080756E+04 # mSQRG_32 
3 3 1.66761356E+06 # niSC!RG_33 



soft SUSY breaking masses 



BLOCK MSU2 Q= 1 . 20195104E+03 # M~2_U 

1 1 1.80692178E+06 # mSURG_ll 

1 2 -5.84935042E-06 # mSURG_12 

1 3 -1.19615520E-02 # mSURG_13 

2 1 -5.84935042E-06 # mSURG_21 
2 2 1.80688177E+06 # mSURG_22 

2 3 4.19949204E+00 # mSURG_23 

3 1 -1.19615520E-02 # mSURG_31 
3 2 4.19949204E+00 # mSURG_32 
3 3 1.23966005E+06 # mSURG_33 



soft SUSY breaking masses 



,OCK MSD2 


Q- 


= 1.20195104E+03 


# M-2_D 


1 


1 


1 


.78917582E+06 


# 


mSDRG. 


.11 


1 


2 


-9, 


.47987536E-06 


# 


mSDRG. 


.12 


1 


3 


5. 


.70618349E-03 


# 


mSDRG. 


.13 


2 


1 


-9, 


.47987536E-06 


# 


mSDRG. 


.21 


2 


2 


1, 


.78914873E+06 


# 


mSDRG. 


.22 


2 


3 


-5, 


.49615942E-01 


# 


mSDRG. 


.23 


3 


1 


5, 


.70618349E-03 


# 


mSDRG. 


.31 


3 


2 


-5 


.49615942E-01 


# 


mSDRG. 


.32 


3 


3 


1, 


.77161275E+06 


# 


mSDRG. 


.33 



soft SUSY breaking masses 



,OCK MSL2 


Q= 


= 1.20195104E+03 


# M~2_L 


1 


1 


2. 


.44436035E+05 


# 


mSLRG. 


.11 


1 


2 


7, 


.11764708E-01 


# 


mSLRG. 


.12 


1 


3 


-1, 


,51306954E+01 


# 


mSLRG. 


.13 


2 


1 


7. 


.11764708E-01 


# 


mSLRG. 


.21 


2 


2 


2. 


.44428213E+05 


# 


mSLRG. 


.22 


2 


3 


7, 


.28699498E+01 


# 


mSLRG. 


.23 


3 


1 


-1. 


.51306954E+01 


# 


mSLRG. 


.31 


3 


2 


7. 


.28699498E+01 


# 


mSLRG. 


.32 


3 


3 


2. 


.41319457E+05 


# 


mSLRG. 


.33 



soft SUSY breaking masses 



,OCK MSE2 


Q= 


= 1.20195104E+03 


# M-S 


!_E 


1 


1 


9, 


,97682891E+04 


# 


mSERG_ 


.11 


1 


2 


6. 


.08801261E-09 


# 


mSERG_ 


.12 


1 


3 


1, 


.01423906E-05 


# 


mSERG_ 


.13 


2 


1 


6, 


,08801261E-09 


# 


mSERG_ 


.21 


2 


2 


9. 


. 97589360E+04 


# 


mSERG_ 


.22 


2 


3 


1. 


. 15378509E-03 


# 


mSERG_ 


.23 


3 


1 


1. 


.01423906E-05 


# 


mSERG_ 


.31 



soft SUSY breaking masses 
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3 2 1.15378509E-03 # mSERG_32 

3 3 9.70150720E+04 # mSERG_33 

# 

BLOCK SuSeFLAVLOWENERGY # PARAMETERS 



1 


3 


.39444376E- 


■03 


# 


Delta rho parameter 


2 


5 


. 11182018E- 


■10 


# 


g_mu - 2 


3 


3 


.73074715E- 


■04 


# 


Br(b -> s gamma) 


4 


1 


.70603079E- 


■13 


# 


Br(tau -> mu gamma) 


5 


7 


.32782350E- 


■15 


# 


Br(tau -> e gsimma) 


6 


9 


.20566951E- 


■17 


# 


BrCm'u -> e gamma) 


7 


5 


.43318338E- 


■16 


# 


Br(tau -> mu m'u mu) 


8 


8 


.82689326E- 


■17 


# 


Br(tau -> e e e) 


9 


6 


.80189300E- 


■19 


# 


Br(mu -> e e e) 



# 

BLOCK FINETUNE # 

1 1.75230800E+02 # delta mZ"2/mZ"2 (mu'2) 

2 8.76436259E+01 # delta mt/mt (mu~2) 

The program prints BLOCK SELMIX in the SLHA output, -which contains the elements of the tree level 6x6 
slepton mixing matrix. The notation and the method of diagonalization is explained in more detail in the 
[Appendix A[ 

7.1.2. Non-SLHA Interface 

The source file for non-SLHA input is runonce . f . To use the traditional SuSeFLAV interface the user should 
execute susef lav by modifying the corresponding input file. This executable takes the following input files, 
sinputs.in for mSUGRA/CMSSM models, sinputs-gmsb. in for GMSB models, s input s-nuhm. in for 
non-universal higgs model (NUHM2) models and sinputs-cnum. in for complete non-universal models. For 
example, to run the program -with complete non-universal boundary conditions type the foUo^wing command 
in the terminal 

./susef lav <sinputs-cnum. in 

The sample input / output below corresponds to 1 . 4 . 4 of [55] with Type-I see-saw mechanism with minimal 
mixing in the context of mSUGRA/CMSSM scenario. To execute the program for the sample point type 
the following command in the terminal 
./susef lav <sinputs.in 



mSUG 


# 


MODEL name 


1 


# 


prinstat (1 for printing output on terminal, otherwise) 


l.OOOOOE+01 


# 


tanbeta 


5.00000E+02 


# 


mO 


O.OOOOOE+00 


# 


aO 


1.05000E+03 


# 


Ml/2 


l.OOOOOE+00 


# 


sgn(mu) 


1.16639E-05 


# 


G_f ermi 


1 . 27934E+02 


# 


alpha_em" (-1) 


1.17200E-01 


# 


alpha_strong 


9.11876E+01 


# 


mz_pole 


2 


# 


one loop or two loops 


1 


# 


quark mixing 


1 


# 


1= rhn on; = rhn off 


CKM 


# 


case: CKM/MNS/USD 


l.OOOOOE-03 


# 


spectrum tolerance 


1 . 72700E+02 


# 


Mtpole 
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4.23000E+00 

1.77000E+00 
l.OOOOOE+06 
l.OOOOOE+09 
l.OOOOOE+14 
O.OOOOOE+00 
O.OOOOOE+00 
O.OOOOOE+00 
O.OOOOOE+00 
O.OOOOOE+00 
O.OOOOOE+OO 
O.OOOOOE+00 
O.OOOOOE+OO 
O.OOOOOE+OO 



# Mbpole 

# Mtaupole 

# MRl 

# MR2 

# MRS 

# Ynu(l,l) 

# Ynu(l,2) 

# Ynu(l,3) 

# Ynu(2,l) 

# Ynu(2,2) 

# Ynu(2,3) 

# Ynu(3,l) 

# Ynu(3,2) 

# Ynu(3,3) 



The Output in SuSeFLAV format is contained in the file susef lav.out. 
******** Begin Program SuSeFLAV ************ 



model 




mSUG 


loop 




2 




tanbeta 




1, 


■OOOOE+Ol 


mO 




5, 


. OOOOE+02 


aO 




0, 


.OOOOE+OO 


M12 




1, 


.0500E+03 


mlO 




5, 


. OOOOE+02 


m20 




5, 


.OOOOE+02 


sign mu 




1, 


. OOOOE+OO 


qmix 




1 




top pole mass 




1, 


.7270E+02 


rhn 




1 




case 




CKM 


MRl 




1, 


.0000E+06 


MR2 




1, 


.OOOOE+09 


MR3 




1, 


.OOOOE+14 



up-type yukawa at high energy: 

2.6401E-06 6.1056E-07 1.2205E-08 

-6.1068E-04 2.6401E-03 1.1119E-04 

3.8672E-03 -1.8624E-02 4.5212E-01 

down-type yukawa at high energy: 

7.9987E-05 -3.1319E-09 7.4082E-08 
-6.2638E-08 1.6000E-03 -7.1355E-06 
3.7936E-05 -1.8270E-04 4.5575E-02 

lepton-type yukawa at high energy: 

1.7539E-05 -1.3061E-10 3.1014E-09 
-2.7282E-08 3.6638E-03 -3.1200E-06 
1.1175E-05 -5.3817E-05 6.4558E-02 



27 



neutrino yukawa at high energy: 



3. 


,1809E- 


■06 


7, 


.3564E- 


■07 


1, 


.4835E- 


■08 


■7. 


,3579E- 


■04 


3, 


. 1810E- 


■03 


1, 


.3396E- 


■04 


4. 


,7106E- 


■03 


-2, 


. 2686E- 


■02 


5, 


.5081E- 


■01 



up-type yukawa at msusy : 



6, 


. 3637E- 


■06 


1, 


.4719E- 


■06 


2, 


.3911E- 


■08 


1, 


.4720E- 


■03 


6 


.3637E- 


■03 


2, 


.6931E- 


■04 


7, 


. 1199E- 


■03 


-3 


.4289E- 


■02 


8, 


.3124E- 


■01 



down-type yukawa at msusy : 



2, 


.4935E- 


■04 


-1, 


.2210E- 


■09 


2, 


.8477E- 


■08 


2, 


.4420E- 


■08 


4, 


. 9870E- 


■03 


-2, 


.7429E- 


■06 


1, 


.4571E- 


■05 


-7, 


.0171E- 


■05 


1, 


.2939E- 


■01 



lepton-type yukawa at msusy : 



2, 


.7626E- 


■05 


-1, 


.0563E- 


■11 


3, 


.7615E- 


■10 


2, 


.2066E- 


■09 


5, 


.7706E- 


■03 


-3, 


. 7838E- 


■07 


1, 


. 3609E- 


■06 


-6, 


. 5535E- 


■06 


9, 


.9328E- 


■02 



Spectrum 



t msusy, q 




1 . 


, ibfvsL+Oo 


alpha_l 




3. 


.6366E-01 


alpha_2 




6, 


.4050E-01 


alpha_3 




1, 


.0297E^H00 


vevl 




2. 


.5253E^H01 


vev2 




2. 


.4230E^H02 


ewtbeta 




9. 


, 5948E^H00 


\mu 




1, 


.2105E^H03 


"gluino 




2. 


.3009E^i-03 


~stop_l 




1. 


, 6702E-I-03 


~stop_2 




1, 


.9889E^H03 


~scharm_R 




2. 


.0550E^i-03 


~scharm_L 




2. 


, 1380E-I-03 


~sup_R 




2, 


,0550E^H03 


~sup_L 




2. 


. 1383EH-03 


sbottom_l 




1. 


, 9606E-I-03 


sbottom_2 




2, 


, 0368E^H03 


s strange _R 




2. 


. 0460EH-03 


sstrange_L 




2. 


, 1394E-f03 


~sd.own_R 




2, 


,0461E^H03 


~sdown_L 




2. 


. 1396EH-03 


~stau_l 




6. 


,3321E+02 


~stau_2 




8, 


.5289E^H02 


~smu_R 




6. 


.4202EH-02 


~smu_L 




8. 


. 6003E-I-02 


~sel_R 




6, 


,4206E^H02 


~sel_L 




8. 


.6012EH-02 


"tausnu 




8. 


.4727E-I-02 


~musnu 




8, 


,5526E^H02 


"elsnu 




8. 


. 5520EH-02 



Higgs Spectrum 
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1.4726E+03 
1.4754E+03 
1 . 1866E+02 
1.4730E+03 
-1.0481E-01 



4.4421E+02 
8.3081E+02 
1 . 2193E+03 
1 . 2280E+03 



8.1458E+02 
1.2179E+03 



mAO 

inh_ charged = 
mhO 
mH 

\alpha_h 

Neutralino spectrum 
Nl 
N2 
N3 
N4 

Chargino spectrum 

CI 
C2 

Low Energy Observables 
Fine tuning, Cmz~2mu~2 
Fine tuning, Cint~2mu~2 
Br(B => s, gamma) 
Br(mu => e, gamma) 
Br(tau => mu , gamma) = 
Br(tau => e, gamma) 
Br(tau => e,e,e) 
Br(tau => mu,mu,mu) 

Br(mu => e,e,e) = 
Br(mu => e in Ti) 
(g_mu - 2) 



flag for the point is AOK 



3, 


.5356E+02 










1, 


. 7680E+02 










3, 


.8430E-04 










3, 


. 1044E-06 


X 


10- 


(■ 


-11) 


5, 


.9654E-06 


X 


10- 


(■ 


-08) 


2, 


.5629E-08 


X 


10- 


(■ 


-07) 


3, 


.0580E-10 


X 


10- 


(■ 


-07) 


1, 


.8321E-09 


X 


10- 


(■ 


-07) 


2, 


.2577E-07 


X 


10- 


(■ 


-12) 


2, 


.0425E-06 


X 


10- 


(■ 


-12) 


1, 


.8491E+00 


X 


10~ 


(• 


-10) 



FLAGS AND THEIR MEANINGS 

AOK = Everything is fine. 

BMUNEG = B_mu is negative at Msusy. 

REWSB = |\mu|~2 < at Msusy. 

MUNOC = Non-convergent |\mu| at Msusy. 

SW2N0C = Non-convergent Sin-2_thetaw at Mz. 

TACSPEC = Spectrum is tachyonic at Msusy. 

NPERTYUK = Non-perturbative yukawa. 

TACSPECMZ = Spectrum is tachyonic at Mz. 

FSNC = Final spectrum non-convergent. 

TACMh = Lightest CP-even neutral higgs tachyonic. 

TACMH = Heaviest CP-even neutral higgs tachyonic. 

TACMA = CP-odd neutral higgs tachyonic. 

TACMHCH = Charged higgs tachyonic. 

VARUNDER = Stepsize is zero while integrating the RGEs 
TACSUP = SUP sector tachyonic. 
TACSDN = SDOWN sector tachyonic. 
TACSLP = SLEPTON sector tachyonic. 
TACSNU = SNEUTRINO sector tachyonic. 
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LEPH = Lightest higgs mass below LEP limit. 
LEPC = Lightest chargino mass < 103.5 GeV. 
LSPSTAU = Lightest stau is LSP . 



8. Outlook 

The MEG experiment has recently reported the latest limits on the rare leptonic decay /i — > e + 7 [5S] . 
Simultaneously, there are results from LHC as well as direct detection dark matter experiments. SuSeFLAV 
is designed to compute supersymmetric spectrum with full flavor violation and further to compute most 
of the leptonic flavor violating observables. The program can be coupled to Dark Matter routines such as 
MicrOMEGAs [37], DarkSUSY [5H] and Superlso [3T] to compute the relic density and direct detection rates. 
The program is free under the GNU Public License and can be downloaded from the following websites: 



http : // cts . iisc . ernet . in/Susef lav/main . html 



http : //pro j acts . hepf orgs . org/ susef lav/ 



In the present version of the program only Type-I seesaw mechanism is implemented. Future versions will 
include other seesaw mechanisms as well as other improvements like inclusion new flavor violating decays 
rates in the quark sector. 
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Appendix A. Tree-Level Masses 

Here we suppress any gauge indices and follow the notation of BPM^^ [70] closely. The Lagrangian 
contains the neutralino mass matrix as 



1 7 
2 



i>°^M-,oi'° + h.c., 



(A.l) 



where ijj" = (B, H^, HuV and 



Ml 






M2 



-Mzcpsw Mzcpcw 



-Mzcpsw Mzspsw \ 

MzCfjCw -Mzspcw 

-/i 

-A* / 



(A.2) 



We use the letter s and c for sine and cosine, so that = sin/3, cp = cos/3 and sw{cw) is the sine (cosine) 
of the weak mixing angle. The 4x4 neutralino mixing matrix is an orthogonal matrix O with real entries. 



Except in our case the sign of fj. parameter is opposite of BPMZ. 
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such that O^Ad^oO is diagonal. The neutralinos Xi are defined such that their absolute masses increase 
with increasing i. Some of theirjiiass values caji_be negative. 

We make the identification = {W^ =F iW^)/V^ for the charged winos and H:^,H^ for the charged 
higgsinos. The Lagrangian contains the chargino mass matrix as 

-lip-^M^J++h.c., (A.3) 
where if+ = {W+, H+f, = (W', H^^Y and 

The chargino masses are found by acting on the matrix A4^+ with a bi-unitary transformation, so that 
U*A4^+V^ is a diagonal matrix containing the two chargino mass eigenvalues, The matrices U and 

V are easily found, as they diagonalizc, respectively, the matrices A^|^_^A^^^ and M.^^_^A4^+. And we have 
taken the U and V matrices to be real. 

The tree-level squark and slepton mass squared values for the family i are found by diagonalizing the 
following mass matrices defined in the {fiL, fiR)^ basis: 



+ ml^ + jsw) ^|c2/3 m„, (^iAu)u - M cot 

TO„, ((ifi)ii - II cot ^) iml)^i + ml^ + ls1yMlc2i) 

{m?Q)u +^1.- {\-\s^) Mlc'2p rrid, {{Ag)u - /xtan/?) 

iT^di ((^d)ii - /utan/3^ (mpii + m| - |s^M|c2^ 

/ {rn\)^i+ml^- {\- s^) Mlc2fi me, ((ig),^ - tan \ 

\ TOe, ((ie)« - A* tan/?) (m|)ji + m^^ - s^M|c2/3 / 



(A.5) 



(A.6) 



(A.7) 



where, w/, e/ are the mass and electric charge of fermion / respectively. The mixing of the first two families 
is suppressed by a small fermion mass, which we approximate to zero. The sfermion mass eigenstates are 
given by 



where Cf is the cosine of the sfermion mixing angle, cos6f, and Sf is sin^j. These angles are given by 



2 m„ I Au — /i cot /3 j 



2 ( — \i tan /3 ) 

tan(20d) = ^ , , ; 1 ^ ^4^72 • (A-10) 



To calculate the lepton flavor violating decays, we diagonalize the full 6x6 sleptonic mass matrix 
{}) by C/^ M] 



(TM?) by C/-^ Al? J7, where U is the sleptonic mixing matrix with real entries. In the gauge basis of 



{sl, fiL,TL,eii, p,ii,fii}'^ the sleptonic mass matrix is defined as 



(AT -II tan P Ij m| + [m^. - s^M|c2^] 1^ 
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where ni|, m| and Ag are the scalar mass matrices and tri- linear coupling matrix respectively. Using the 
full 6x6 sleptonic mixing matrix we define the lepton flavor violating couplings following Hisano et al. |14j . 
With these couplings we calculate the rare lepton flavor violating decays as described in section [6| 

Given values for tan/3 one can write the CP-odd Higgs-boson mass, m^i, and the other Higgs masses at 
tree level by 

2B., . o ^^^^2) 

'H.h - o I '"-A ^ V V'-A ^ "^Zi ~ '±"'-A^"Z^2S I ' (A.13) 

and 

m]j± = m\ + . (A. 14) 

The CP-even gauge eigenstates {H^, H^) are rotated by the angle a into the mass eigenstates (iJ, h) as 
follows, 

' Co 





(A.15) 



At tree level, the angle a is given by 



M 



2 



tan2a = ""^ ' ^' f tan 2/3 . (A.16) 



Appendix B. One Loop Threshold Corrections 

We compute flavor conserving complete one-loop corrections to masses and couplings in MSSM following 
BPMZ [7(3J. In this appendix we present the summary of one-loop corrections implemented in SuSeFLAV. 

• Quarks and Leptons 

The corrections to fermions are evaluated at the weak scale. The running masses are related to 
the corresponding DR masses by the expression 

= m^.-Sf^*^^(TO2) (B.l) 

Where, 'S^^^^^ {nij) is the one loop self-energy of the fermion /. We follow equation D.18 of BPMZ 
[70] . At the present version of SuSeFLAV we add the correction only to the third generation fermions. 

• W and Z Bosons 

Corrections to W and Z bosons are evaluated at the weak scale as well as EWSB scale. The DR 
running W and Z boson mass squared at the scale Q (Mf (Q) and Myy{Q)) are related to the physical 
pole mass of gauge bosons as follows 

M| = MUQ) - ni^(M|) (B.2) 
= MUQ) - ^wwiMlr) (B.3) 

Where, H^^{M^) and n^^(M^) are the transverse part of self energy terms. Consult appendix D 
of BPMZ 70 for detailed discussion and the complete expression of the self energy terms. 

• sQuarks and sLeptons 
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Flavor conserving one-loop masses and mixings for squarks and sleptons are evaluated by diagonalizing 
M^^{p^) for eigenvalues at the EWSB scale. 



M 



fnj' fRfrV-f / 



JLJR 
JRJR 



n 



iti 



M?^ is the 2x2 tree level mass matrix defined in equations 



of the above matrix 



ix ( n ? r 

\ JLIL 



^/fli-R. 



A.5 



A.6 



RiP') 
^iRfRiP') 

and 



(B.4) 



A.7 



The self-energy terms 



are defined in appendix D of BPMZ 



• Higgs Bosons 

The mass of the two loop CP-even higges are obtained by diagonalizing A4^{p^) at the EWSB scale. 
Where the matrix A4^{p^) is given by, 



Mlcl + m\sl - n,^,^ {p^) + ti/vi -(M| + rri\)spcp - H,^,, {p' 

f2 I -2 



,(p')+t2/w2 



(B.5) 



The self energy terms n^^sj, IIs^s^, IIs^s^ and H^^s^ are presented in appendix D of BPMZ [70] . 
Whereas, one loop tadpole contributions ti/vi and ^2/^2 are presented in appendix E of BPMZ [70] . 
Besides the complete one loop correction to m/j and mjj , we evaluate the top mass enhanced dominant 
two loop corrections provided in |71j . From version 1.2 onwards, we have implemented the full two 
loop Higgs corrections due to Slavich et. al [72] . 

The tree level mass of CP-odd higgs boson, itla and charged higgs boson mu± is given by equations 
|A.12| and A. 14 respectively. The one loop correction to masses is evaluated at EWSB scale and given 

by 



rhH± 



(B.6) 
(B.7) 



Where, 6^ = s} hhi + c| ^2/^2- 
• Neutralino and Chargino 

The one loop corrected neutralino mass matrix has the following form. 



+ \ (5M^o{p') + 5MlAp')) 



where 



^Up' 



(B. 



(B.9) 



Where, M^o is the tree level neutralino mass matrix defined in A. 2 And the matrix 6A4^o{p ) is the 
one loop correction to the neutralino mass matrix. 

The one-loop chargino mass matrix is as follows. 



(B.IO) 



Where, ^A^+ is the tree level chargino mass matrix defined in 



A.4 



We evaluate the self energies 



^LR s with p = (3(EWSB scale). See appendix D of BPMZ for the complete expressions of these 
self-energies. One loop neutralino and chargino mass matrices are then diagonalized to obtain the 
eigenvalues and eigenvectors at EWSB scale. 
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• Gauge Couplings and sin^ 9w 

The correction to DR electromagnetic coupling at the weak scale is given by, 

= — T-r and aem = ^^^tttttt (B.ll) 

l-Aa 127.934 ^ ' 

Where Ad contains contribution from SM particles as well as SUSY particles. Similarly, the correction 
to strong coupling as at the weak scale is given by, 

(B.12) 



1 - Aa, 

Aofs receives contribution from top quark, gluino and squarks. The exact expression given by equation 
3 of BPMZ [70] . The correction to DR weak mixing angle is given by, 

A'jrry 

sin^ 29w = r- n (B.13) 

y2M|G^(l-Af) ^ ' 

The exact expression of Ad and Af is presented in the appendix C of BPMZ [70] . 

• Gluino 

One loop corrected gluino mass uig at scale Q is given by, 

rrig = AhiQ) - I]g(mg) (B.14) 

where A'hlQ) is the tree level gluino mass generated by the RGEs and Tig{mg) is the one loop self- 
energy contribution to the gluino mass. See equation D.44 of BPMZ [70] for the complete expression. 



Appendix C. Program File Structure 

In this appendix we briefly list the subroutines of interest in the program and its purpose (for more 



details see the technical manual provided with the program). Fig. (C.3| depicts the program structure and 
the flow of variables in the program. The subroutine name is printed first followed by the respective file 
which contains the subroutine. 

The program begins by reading inputs from the main program contained in the files 'runonce.f, 
'runslha.f ' (for SLHA input format), 'scEinning.f '. The convention followed in the program is such that all 
the input parameters which do not receive one-loop susy threshold correction such as the masses of the first 
two generation of quarks and leptons are contained in the file 'stdinputs .h'. Whereas, parameters such as 
sm^Owi gauge couplings and the third generation quark and lepton masses which receive one- loop SUSY 
threshold correction are contained in the file 'SuSemain.f '. Moreover, the complete 3x3 structure of the 
Yukawas and soft terms are defined in 'SuSemain.f. 



• rgeiterate . f : The heart of the program. The iterative algorithm is contained here. 

• sof tspectrum. f : Given the RGE output at low scale, this routine computes the complete tree level 
SUSY spectrum by diagonalizing the mass matrices. 

• slha.f : Contains the ingredients for SLHA input/output interface. 

• mueconvernew. f : Contains routines which calculate the decay rates and branching fractions of rare 
lepton flavor violating decays. 

• mssmrge.f : The complete MSSM two loop RGEs are contained in this file. 

• oB.e\oo^particle . f : A set of several files which contain the one loop corrections to the corresponding 
particle^ which could be gauge, neutralino, charging, sfermion etc. 



34 





riuiouce 




runonce.f 


Program • 


ruiislha 


runslha.f 




scaiispace 




scanning.f 



Subroutine 
coinpletetruu 
spectrumtl . f 



Subroutine suseflav 
SuSemain . f 



Recursive Sub- 
routine rgeit 
rgeiterate . f 



Subroutine 
print slha 
slha.f 



Subroutine 
bsg bsg . f 
Subroutine 
muecorvernew 
mueconvernew . f 



Subroutine nissnu'un 
runrges . f 



Subroutine 
softspectruni 
softspectrnm. f 



Recursive Subroutine 
iterate ewsbiterate . f 



Subroutine tad- 
polel, tadpole2 
onelooptadpole . f 



Subroutine rewsbcor 
ewsbiterate . f 




oneloopself energy . f oneloophiggs . f 

oneloopsf ermion . f oneloopchargino . f 
oneloopneutralino . f twoloophiggs . f 



Subroutine runtomz 
spectrumtl . f 



Subroutine nissni_niz 
spectrumtl . f 



Subroutine 
softspectrnm 
softspectrnm. f 



Subroutine topcor, 
bottomcor. taucor 
oneloopf ermion. f 



Figure C.3: File structure and program flow 



• oneloopPV.f : The analytical form of all scalar one loop Passarino-Veltman functions are contained 
in this file. 

• math . f : Routines for matrix manipulations, integration routine for ordinary differential equations 
using Runge-Kutta with adaptive step size and other general purpose functions such as random number 
generator. 

Appendix D. Spectrum Comparison with Other Codes 

In this section we compare the spectrum of SuSeFLAV with other codes for two different parameter space 
points form [SB]. For both the points we have set RHN equals to zero. 
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Table D.2: Spectrum comparison for point 10.1.5 of [96| 



10.1.5 


mo 


= 175, Ml ^ 700, tan/3 

2 


= 10, Aq = 0, sign(^) > 


Parameters 


SuSeFLAV 




SOFTSUSY 


SuSpect 




SPheno 


W+ 


8.05096462E^ 


-01 


8.03860774e-t-01 


8.04813477EH 


hOl 


8.03115948E+01 


h 


1.16294107E-( 


-02 


1.16694326e-|-02 


1.16799828EH 


h02 


1.16829168E+02 


H 


9.66133037E4 


-02 


9.74430227e-h02 


9.68302501EH 


h02 


9.79779724E+02 


A 


9.65805171E^ 


-02 


9.74206243e-t-02 


9.68101504EH 


h02 


9.79577186E-H02 


H+ 


9.69592528E^ 


-02 


9.77728714e-t-02 


9.71660981EH 


h02 


9.83266787E-t-02 


QL 


1.43474627E^ 


-03 


1.43580606e-H03 


1.43214020EH 


h03 


1.43992830E-t-03 


QR 


1.37749634E^ 


-03 


1.37825389e-t-03 


1.37724979E-I-03 


1.38138749E+03 


h 


1.11659996E^ 


-03 


1.11105612e-|-03 


1.11462082EH 


h03 


1.10949455E-I-03 


h 


1.35668185E^ 


-03 


1.35225046e-|-03 


1.35767238EH 


h03 


1.35426026E-I-03 


hi 


1.31765879E^ 


-03 


1.31707770e-t-03 


1.32640786EH 


h03 


1.31861062E-t-03 




1.36683066E^ 


-03 


1.36729669e-l-03 


1.36701474EH 


h03 


1.37028335E-)-03 


II 


5.01303257E^ 


-02 


5.01757612e-H02 


4.94643201EH 


h02 


5.01101877E-I-02 


Ir 


3.23056643E+02 


3.18458635e-t-02 


3.13310890E-I-02 


3.18342964E-t-02 


n 


3.15781781E+02 


3.11269901e-K02 


3.06660771E-I-02 


3.11160940E-I-02 


T2 


5.01464374E^ 


-02 


5.01745556e-t-02 


4.94979107EH 


h02 


5.01126010E-I-02 


i>i 


4.94057805E^ 


-02 


4.95324112e-H02 


4.88544346EH 


1-02 


4.94578335E-f02 


i>T 


4.92669260E^ 


-02 


4.93751173e-H02 


4.87265068EH 


h02 


4.92997226E-t-02 


9 


1.56391437E^ 


-03 


1.56476030e-t-03 


1.56495608EH 


h03 


1.57232851E-H03 


X? 


2.95219584E+02 


2.91289971e-f-02 


2.92926561EH 


h02 


2.92509855E-I-02 


x§ 


5.45858411E^ 


-02 


5.51079728e-h02 


5.52119478EH 


h02 


5.52386977E-t-02 


xt 


8.44975403E^ 


-02 


8.50445968e-h02 


8.60389258EH 


h02 


8.56133860E-I-02 


xl 


8.56487466E4 


-02 


8.62039492e-|-02 


8.72046145EH 


h02 


8.67579765E+02 


xf 


5.35344839E+02 


5.51280182e-l-02 


5.52043599E-I-02 


5.52530182E-H02 


xt 


8.49067665E+02 


8.61668017e-K02 


8.71711677EH 


h02 


8.67856046E-I-02 




8.37298648E^ 


-02 


8.44864881e-l-02 


8.55544693EH 


h02 


8.50798325E-1-02 


-^SUSY 


1.19242602E^ 


-03 


1.19031414e-t-03 


1.19174751EH 


h03 


1.22578828E-H03 


Mgut 


1.26703749E^ 


-16 


1.67323027e-M6 






1.53713886E+16 
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Table D.3: Spectrum comparison for point 40.2.3 of I96| 



40 . 2 . 3 


mo = 


650, Ml = 550, tan/3 


^ 40, ^0 -500, sign{^i) > 




Parameters 


SuSeFLAV 




SOFTSUSY 


SuSpect 




SPheno 




W+ 


8.05066107E^ 


-01 


8.03827802e-|-01 


8.04773725EH 


hOl 


8.03087734E-t-01 


h 


1.16958052E^ 


-02 


1.17297317e-h02 


1.17322254EH 


h02 


1.17682268E^ 


h02 


H 


7.30962046E4 


-02 


7.47689856e-H02 


7.18511292EH 


h02 


7.72497154E- 


h02 


A 


7.31051866E4 


-02 


7.47643928e-l-02 


7.18551942E-I-02 


7.72499188E^ 


h02 


H+ 


7.35889065E+02 


7.52307728e-H02 


7.23362920E-H02 


7.77278228E-I-02 


UL 


1.30348659E^ 


-03 


1.28439273e-t-03 


1.30388612EH 


h03 


1.31020706E-I-03 




1.26813307E^ 


-03 


1.24781930e-H03 


1.27013696EH 


h03 


1.27389691E- 


1-03 


dL 


1.30579237E^ 


-03 


1.28668811e-l-03 


1.30621620E-I-03 


1.31246645E- 


h03 


dR 


1.26498184E^ 


-03 


1.24450694e-|-03 


1.26672885EH 


h03 


1.27054542E^ 


h03 


ii 


9.24983971E^ 


-02 


9.21498284e-|-02 


9.69783691EH 


h02 


9.19586596E^ 


h02 


h 


1.14470194E^ 


-03 


1.14136271e-h03 


1.17704096EH 


h03 


1.14412090E-. 


h03 


h 


1.08597531E^ 


-03 


1.08516841e-l-03 


1.11715522E-I-03 


1.09021482E-I-03 


h 


1.15972419E^ 


-03 


1.16046164e-H03 


1.19253887EH 


h03 


1.16721088E- 


h03 


II 


7.46153736E^ 


-02 


7.46218781e-t-02 


7.43808878EH 


h02 


7.46077658E- 


h02 


Ir 


6.83925820E^ 


-02 


6.82393714e-}-02 


6.80821036EH 


h02 


6.82370617E- 


1-02 


n 


5.23418984E^ 


-02 


5.18531351e-t-02 


5.54473085EH 


h02 


5.18371511E- 


h02 




6.95666250E^ 


-02 


6.94111081e-H02 


7.16257580EH 


h02 


6.94363534E^ 


h02 


I'l 


7.40378630E^ 


-02 


7.41675551e-H02 


7.39679081EH 


h02 


7.41481196E^ 


h02 




6.79971865E^ 


-02 


6.77798789e-t-02 


6.94510135EH 


h02 


6.77596601E-t-02 


9 


1.27460200E+03 


1.27746572e-|-03 


1.28080280EH 


h03 


1.28539752E^ 


h03 


x\ 


2.33065416E^ 


-02 


2.29995196e-|-02 


2.31351813EH 


1-02 


2.30719919E^ 


1-02 


xl 


4.32028446E^ 


-02 


4.38269965e-f-02 


4.43044390EH 


h02 


4.39375323E^ 


h02 


xl 


7.36025983E^ 


-02 


7.45547839e-t-02 


9.13965835EH 


h02 


7.54411085E- 


h02 


xl 


7.44572370E+02 


7.54116707e-t-02 


9.19380416E-I-02 


7.62689372E-t-02 


xf 


4.22498882E+02 


4.38347771e-t-02 


4.43031333EH 


h02 


4.39557714E^ 


h02 


xt 


7.40879408E^ 


-02 


7.54555638e-|-02 


9.20007967EH 


h02 


7.63794652E^ 


h02 




7.29597776E^ 


-02 


7.41004045e-f02 


9.12311446EH 


h02 


7.50392113E^ 


h02 


-^SUSY 


9.97283261E^ 


-02 


9.96175665e-H02 


1.03350532EH 


h03 


1.02572777E- 


h03 


Mgut 


1.23171917E4 


-16 


1.88768830e-t-16 






1.75795520E^ 


hl6 
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